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Chapter  1 

INTRODUCTION 


The  very  simplicity  of  the  Hall  thruster  structure  with  no  screens  or  other  control  surfaces, 
makes  it  difficult  to  understand  the  interrelationships  which,  in  the  end,  localize  and  shape 
the  various  plasma  regions  existing  in  the  accelerating  channel.  Bishaev  and  Kim[l]  distin¬ 
guish,  in  their  experiments  i)  an  upstream  diffusion  region,  with  negligible  electric  field,  no 
ion  production,  predominant  ion  motion  towards  the  anode,  and  low  plasma  temperature;  ii) 
an  intense  ionization  layer,  placed  near  the  channel  exit,  with  peaks  of  plasma  density  and 
temperature;  and  iii)  an  acceleration  region,  where  most  of  the  potential  drop  takes  place,  and 
extending  outside  the  channel  into  iv)  the  plume.  These  characteristics  of  the  plasma  flow 
are  confirmed  by  other  experimental  observations [2,  3]  and  by  two-dimensional(2-D)  plasma 
simulations[4].  Since  the  radial  magnetic  field  is  usually  shaped  with  a  peak  near  the  channel 
exit,  the  above  plasma  structure  has  often  been  explained  as  simply  a  reflection  of  the  magnetic 
field  distribution. 

V.  Kim,  in  a  recent  review  paper  on  Hall  thruster  physics[5],  recognizes  that  there  is  still 
no  model  describing  adequately  the  plasma  dynamics  inside  the  accelerating  channel.  The 
importance  of  developing  a  reliable  model  that  improves  our  understanding  of  the  factors  which 
control  plasma  structure  becomes  evident  when  we  attempt  to  improve  upon  existing  Hall 
thruster  designs,  particularly  regarding  (i)  thrust  efficiency,  (ii)  erosion,  and  (iii)  control  of 
oscillations.  Present  design  criteria  are  based  mainly  on  testing  experience.  Erosion  of  the 
thruster  walls  is  the  main  factor  limiting  the  thruster  lifetime.  The  erosion  rate  is  governed 
by  local  plasma  density,  electron  temperature,  and  sheath  drop,  and  each  of  these  may  show 
order-of-magnitude  variations  within  the  channel  making  it  crucial  to  have  theoretical  guidance 
for  steering  their  peaks  to  the  least  damaging  locations  possible.  Plasma  self-oscillations  are 
manifested  by  a  modulation  of  the  discharge  current,  which  sometimes  is  fairly  deep,  but  at 
other  times  it  is  much  weaker  or  almost  absent,  and  the  reasons  for  these  widely  different  and 
badly-predictable  behaviors  are  not  understood  at  present. 

One-dimensional(lD)  macroscopic  models  can  be  simple  but  powerful  tools  to  understand 
the  physics  and  parameters  that  govern  the  thruster  response  and  to  determine  the  basic  scaling 
laws  for  thruster  magnitudes.  Furthermore,  these  fast-to-solve  models  can  be  the  ideal  comple¬ 
ment  to  time-consuming  PIC-MCC  2D  simulations[6,  7],  for  Hall  thruster  design  and  analysis. 
In  general,  published  ID  analyses  (Refs.  [8,  9,  10,  11]  are  some  of  them)  put  the  emphasis  on 
achieving  satisfactory  performance  prediction  which,  if  a  few  reasonable  assumptions  are  intro¬ 
duced  about  the  flow  structure,  is  not  too  difficult  to  achieve  (particularly  regarding  thrust). 
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However,  these  models  fail  to  describe  the  whole  plasma  structure  described  by  Bishaev  and 
Kim,  and,  in  general,  the  model  formulation  was  not  consistent  enough. 

In  1998,  E.  Ahedo  and  M.  Martinez- Sanchez  decided  to  develop  and  analyze  a  macroscopic 
model  gathering  reliability  and  clarity.  Quantitative  accuracy  in  the  predictions  should  not  be 
the  main  goal.  Instead  emphasis  will  be  given  to  the  consistency  of  the  governing  equations 
and  their  boundary  conditions,  and  to  include  enough  basic  physics  to  illuminate  the  complex 
interplay  of  phenomena  which  shape  the  Hall  thruster  plasma.  Preliminary  results[12,  13]  were 
promising  and  constitute  the  seed  of  the  present  Project. 

We  have  organized  the  main  body  of  this  work  in  six  chapters.  In  Chap.  2  we  discuss  the 
general  formulation  of  the  model,  from  the  discussion  of  the  phenomena  to  be  included  in  the 
plasma  equations  and  their  function  form  to  the  geometry  and  unavoidable  simplifications  (in 
the  plume  treatment,  for  instance).  In  Chap.  3  we  discuss  the  basic  model,  that  is  the  model 
we  think  ’minimum’  to  reproduce  the  rich  plasma  structure  described  by  Bishaev  and  Kim. 
In  Chap.  4  we  discuss  a  supplementary  model  of  the  radial  interaction  of  the  plasma  with 
the  lateral  walls  of  the  thruster.  In  Chap.  5  we  analyze  the  effects  of  the  radial  plasma/wall 
interaction  on  the  basic  steady-state  solution  of  Chap.  3.  In  Chap.  6  we  solve  a  model  that 
includes  heat  conduction  effects  (but  no  radial  effects).  Finally,  In  Chap.  7  we  approach  the 
problem  of  self-oscillations  from  a  stability  analysis  of  the  steady-state  plasma  structure  of 
Chap.  5. 


Chapter  2 

FORMULATION  OF  THE  1-D  MODEL 


2.1  Model  hypotheses  and  geometry 

Geometrical  sketches  of  the  thruster  and  the  main  features  of  the  model  are  drawn  in  Figs. 
2.1(a)  and  2.1(b).  The  macroscopic  model  considers  the  plasma  to  be  composed  by  three  inde¬ 
pendent  fluids:  electrons  (e),  ions  (i),  and  neutrals  (n),  and  formulates  equations  for  particle, 
momentum,  and  energy  conservation  of  each  species. 

The  channel  is  cylindrical  of  length  L  and  has  a  constant  radial  area  Ac  (and  radial  thickness 
he);  if  X  is  the  axial  variable,  we  take  xa  =  0  at  the  anode,  and  xe  =  L  at  the  channel  exit.  The 
plume  is  modeled  as  a  divergent  cylindrical  jet  with  a  radial  area  A{x)^  which  will  be  determined 
as  part  of  the  solution.  Electrons  are  injected  into  the  plasma  plume  through  a  virtual  cathode 
or  neutralization  surface  (point  P)  placed  at  a  distance  Lep{=  xp  —  xe)  from  the  channel  exit. 
The  voltage  difference  between  points  A  and  P  is  the  discharge  voltage  Va  =  (j)A  —  (j)p- 
The  applied  magnetic  field  is  assumed  radial 

B{x)  =  B{x)er,  (2.1) 

and  its  strength  is  such  that  ions  can  be  assumed  unmagnetized  whereas  electrons  are  strongly 
magnetized,  and  follow  a  quasi-closed  drift  motion  across  magnetic  field  lines. 

The  electron  current  delivered  at  the  neutralization  surface  is  the  discharge  current  Id- 
One  part  of  this  current  diffuses  inwards  (into  the  near-plume  and  the  channel)  and  ionizes 
the  mass  flow  of  neutrals,  m,  injected  at  the  anode.  The  other  part  of  the  electron  current 
flows  outwards  (into  the  far  plume)  and  neutralizes  the  ion  current;  subscript  oo  will  refer  to 
downstream  conditions  far  away  from  the  cathode.  To  facilitate  this  flow  of  electrons  out  of  the 
cathode  the  voltage  profile,  ^(x),  is  expected  to  be  minimum  at  point  P.  The  voltage  difference 
4^00  —  (f^p  >  0  can  be  interpreted  as  the  cost  of  the  electron  emission. 

Since  the  Debye  length  is  generally  very  small  (of  the  order  of  10  —  100//m)  compared  to 
the  channel  length,  the  plasma  is  considered  quasineutral  everywhere  except  in  an  electron- 
repelling  sheath  attached  to  the  anode;  this  is  region  AB  in  Fig.  2.1(b),  with  point  B  marking 
the  transition  to  the  sheath.  In  this  sheath,  a  self-adjusted  potential  bias, 

4>sh  =  —  <pA  >  0., 

equates  the  electron  thermal  flow  collected  at  the  anode  to  the  quasineutral  diffusion  flow  of 
electrons  across  the  thruster  channel.  We  consider  that  only  neutral  gas  is  supplied  externally 
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through  the  anode.  Then,  since  the  anode  sheath  is  thin  and  ion-attracting,  ions  must  flow 
into  the  sheath  from  the  quasineutral  region  adjacent  to  the  anode.  The  presence  of  a  region 
with  reverse  ion  flow  between  the  anode  and  the  ionization  layer  is  well  reported  in  Refs.  [1] 
and  [5].  The  reverse  ion  flow  is  recombined  at  the  anode,  and  re-emitted  as  neutral  gas  to  the 
channel,  adding  itself  to  the  gas  flow  supplied  externally. 

The  interaction  of  the  plasma  discharge  with  the  lateral  walls  of  the  thruster  makes  the 
plasma  flow  two-dimensional  (2D).  Instead  of  tackling  directly  this  2D  problem,  which  would 
require  a  sophisticated  numerical  procedure,  our  proposal  is  to  solve  separately  two  approxi¬ 
mate  ID  problems,  for  the  axial  and  radial  directions,  and  then,  to  couple  them  with  adequate 
parameters  and  variables.  The  main  ID  axial  model  will  use  variables  averaged  radially,  and  in¬ 
cludes  the  interaction  of  the  plasma  with  lateral  walls  through  source  terms  for  production/loss 
of  particles,  momentum  and  energy.  The  ID  radial  model  must  provide  the  way  to  compute 
these  source  terms.  Of  course,  this  ID  +  ID  procedure  is  not  exact  but  presents  the  advantages 
of  (i)  to  reduce  the  problem  to  conservation  and  ordinary  differential  equations,  (ii)  to  lead  to 
a  quicker  resolution,  and  (iii)  to  provide  a  lot  of  physical  insight. 


2.2  Plasma  macroscopic  equations 


From  the  Boltzmann  equation,  the  three  basic  macroscopic  equations  for  a  given  plasma  species 
(species  subscript  is  omitted)  are: 

i)  Particle  conservation, 

^  +  V-nv  =  S.  (2.2) 

C/6 

ii)  Momentum  conservation, 

Otiv 

m— - h  mV  •  (nvv)  =  —  Vp  +  ZeniE  +  v  A  B)  +  M.  (2.3) 

ot 

iii)  Energy  conservation. 


9/1  2 


■nT^  +  V  •  (^^nmv^  +  ^nT^v  -h  V  •  g  = 


-  V  ■  pv  +  Zenv 


where:  v  is  the  macroscopic  velocity;  q  is  the  conductive  heat  flux;  S,  M,  and  Q  represent 
external  sources  and  collisional  processes  with  other  species;  and  other  symbols  are  conventional. 
To  complete  these  macroscopic  equations,  a  transport  equation  for  q,  and  expressions  for  the 
three  source  terms  must  be  specified. 

From  the  above  fundamental  equations  related  equations  of  interest  are: 

-  Velocity  evolution. 


f  —  -h  t;  •  v)  V  =  —  Vp  -h  Zen(E  -\-  v  A  B)  -\-  M\ 
\dt  ) 


mn 


with  M'  =  M  —  Smv. 

-  Mechanical  energy  evolution, 

d  1 


n 


f  —  -h  v  •  v)  -mv^  =  —V  ■  Vp  -h  Zenv  ■  E  +  v  ■  M  —  Smv^. 
\dt  /  2 


(2.5) 


(2.6) 
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-  Mechanical  energy  conservation, 

f  ^nmv^)  +  V  •  f  ^nmv^v)  =  —v  •  Vp  +  Zenv  •  E  +  v  •  M  —  ^Smv^. 

-  Internal  energy  conservation, 

^(^nT)  +V-  +  V  •  q  = -pV  •  r;  +  Q', 

with  Q'  =  Q  —  V  ■  M.  T  ^Smv^. 

Zi 

-  Internal  energy  evolution. 


n 


(2.7) 


(2,8) 


(2.9) 


The  set  of  conservation  equations  for  all  plasma  species  is  completed  with  Poisson  equation 
for  the  electric  field  E  or  the  electric  potential  (j), 


cqV  ■  E  =  ZaTia,  with  E  =  —V<j). 


(2.10) 


2.2.1  The  2D  axisymmetric  equations 

Assuming  azimuthal  symmetry,  d/dO  =  0,  and  taking  into  account  that  the  magnetic  field 
is  radial,  the  equations  for  the  density,  the  velocity  and  the  internal  energy  become 


dn  d  Id  „ 

—  +  —nv^  +  -^nrvr  =  A, 
ot  ox  r  or 


d 

mn  I  —  +  v^—  +  ■Wr—  I  Vx  =  - 


'_a 

.dt 

d 


'  dx 
d 


d  \ 
^"dr) 


dnT 


drJ~^  dx 
d 


+  Zen{E^  -  Bve)  +  M^, 


(U  U  0  V  \ 

+  v^—  +  Vr^  +  —  ]ve  =  ZenBv^  +  Mg, 
dt  dx  dr  r  / 

( d  d  d  \  Vg  dnT  „  „  , ,, 

+  '^r^)vr - =  — ^ - h  ZenEr  +M^, 

\dt  dx  dr/  r 


d  /3 


d  /5 


dr 
1  d  .5 


d 


Vr  d 


dt^2 


r  dr  ^2 


dx 


r  dr 


(2.11) 

(2.12) 

(2.13) 

(2.14) 

(2.15) 


2.2.2  The  approximate  ID  axial  model 


This  is  defined  by  assuming  plasma  variables  independent  of  r  and  including  the  effects  of 
radial  boundaries  as  extra  source  terms  in  the  axial  equations.  From  Eqs.  (2. 11)- (2. 15),  the 
axial  equations  for  the  radially-averaged  magnitudes  (represented  with  a  hat)  are 


dh  a  _  _  Q  ,  Q 
"57  +  —  S  +  S^, 

dt  dx 


mn 


{  o  ^  o  \  ^ 

la = 


dnT 

dx 


+  Zen{Ex  —  Bvg)  +  M'^  -\-  M'^^, 


M 

d  . 
mn  I  —  -h 
.dt 


( dx 


]ve 


ZenBv^  +  Mg  +  M'^g, 


(2.16) 

(2.17) 

(2.18) 
(2.19) 
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where  S^,  M'^^^  ^'wei  Q'w  represent  the  wall  source  terms.  Notice  that  Eq.(2.14)  is  not  used  in 
the  axial  model. 

In  the  cylindrical  thruster,  the  convenient  radial  average  for  a  generic  magnitude  /  seems 
to  be 

f{x,t)  =<  rf{r,x,t)  >=  2  ^  2  /  (2.20) 

where  r\  and  r2  are  the  inner  and  outer  radius,  respectively.  Then,  the  wall-source  terms  take 
the  form 


d 

<  —rnvr  >, 
or 


-  d 

M-ujx  =  -  <  -^{rmnVrVx)  >, 

or 


\r  (jr-^ 


(2.21) 

(2.22) 

(2.23) 

(2.24) 


2.3  Collisional  processes 


We  discuss  here  the  different  collisional  phenomena  among  the  plasma  species  and  the  contri¬ 
butions  they  yield  to  terms  S,  M,  and  Q',  in  Eqs.(2.16)-(2.19).  Wall-source  terms,  Eqs.(2.21)- 

(2.24)  will  be  discussed  in  Chap.  4. 


2.3.1  The  ionization  rate 

We  consider  that  (i)  ionization  is  produced  by  electron-neutral  collisions  exclusively,  and  (ii) 
we  neglect  double  ionization  (i.e.,  Zi  =  1).  Here  we  look  for  an  analytical  expression  of  the 
ionization  source  term,  which  reproduces  correctly  enough  its  dependence  with  the  electron 
temperature. 

The  ion  production  can  be  expressed  as 

Si  =  Uel^i  l^i  =  n„i?j(Te), 

where  Ri{Te)  is  the  effective  ionization  rate.  Since  the  dependence  of  Ri  on  the  electron  tem¬ 
perature  is  going  to  be  crucial  in  shaping  the  ionization  region,  it  is  important  to  use  a  correct 
expression  for  it. 

Assuming  that  the  electron  motion  satisfies 

\ve\  < 


the  electron  distribution  function  can  be  taken  local-Maxwellian, 


(2.25) 
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where  E  is  the  electron  energy.  Then,  the  ionization  rate  can  be  written  as 

o  poo 

Ri{T,)  =  —  /  ai{E)Ef,{E)dE, 

Jo 

where  ai{E)  is  the  ionization  cross-section. 

To  work  with  simple  analytical  expressions,  ai{E)  (for  xenon)  is  given  by 

ai{E)  =  (Tjo  1  -  exp{E/Ei  -  1) 

Ei  =  12.1  eV,  (Tio  =  5  X  10 


-20  ^^2 


(2.26) 


(2.27) 


for  E  >  Ei,  and  ai{E)  =  0  for  E'  <  Ef,  E^  is  the  energy  for  primary  ionization.  Figure  2.2(a) 
compares  this  expression  with  experimental  values,  showing  that  it  is  a  good  approximation  for 
electron  energies  E  up  to  lOOeV,  which  is  the  range  of  interest  for  a  Hall  thruster. 
Substituting  fe{E)  and  ai{E)  in  Eq.  (2.26),  one  obtains 

Ri{Te)  =  di{Te)ce,  with  dj(Te)  =  aio(l  +  (  ~ 

and  Ce  =  ^/STe/Tvme.  Figure  2.2(b)  compares  this  expression  of  di(Te)  with  the  one  obtained 
numerically  from  experimental  values  of  ai{E).  The  agreement  is  very  good  up  to  Tg  ~  40eV 
and  satisfactory  up  to  Tg  ~  100  eV.  Although  it  is  easy  to  obtain  better  improvements  with 
more  accurate  expressions  of  (Ti{E),  we  find  it  worthless  for  the  objectives  of  the  present  model. 
The  ionization  frequency  can  be  computed  from 


- rin - ^  ^  - - —  ^  ^  ® 


1020  j^-3  Vl  eV^'  5  X  10“^°m2 
The  neutral  density  of  xenon  (m^  ~  2.19  •  10“^^kg)  follows 

m  10  cm^  300  m/s 


m 

n„.  - - 


X 


rriiAvn  1  mg/s  A 
For  the  mass  flow,  the  equivalences 


X 


X  1.47  X  10^®  m-^ 


131  gr  cm^ 

1  mg/s  =  0.731  amp-eq,  1  seem  ~ - - - ^  ~  0.0975  mg/s, 

^  22.4  ht  mm  ’ 

can  be  useful. 

Ionization  produces  a  momentum  exchange  between  ions  and  neutrals  given  by 


(2.29) 


(2.30) 


2.3.2  Charge- exchange  collisions 


Charge-exchange  collisions  between  ions  and  neutrals  produce  a  momentum  exchange  given  by 

AL  gx  —  '^n)j 
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where  the  charge-exchange  source  term  can  be  expressed  as 


Sex  ^cx  '^nCn^cxi 


(2.32) 


with  Cin  =  \vi  -  Vn\  and  Ucx  given  by[14] 


a  ex  —  1  —  0.1  In 


.  1  km/sv 


X  8.16  X  10“^W. 


(2.33) 


Computing  the  frequency  for  charge- exchange  from 

Cjr*  ^cx  r  1  n4  — 1 

^cx  =-  ?  X  - - -  X  - -  X  5  X  10  s  , 

1020  ikm/s  5  X  10“^®m2 

one  finds  that  it  is  a  secondary  effect  which  can  be  neglected  in  the  fundamental  model  of  the 
plasma  flow. 

2.3.3  Electron-neutral  elastic  collisions 

The  electron-neutral  collision  frequency  is  expressed  as 

^en 

with  aen  depending  on  Tg.  Since  these  collisions  are  important  only  in  the  upstream  part  of 
the  channel  we  just  take 

aen  —  const  =  27  X  10“^°  m^. 

Using  this  average  value  one  has 


7/  _  V 

-  1020  m-3 


Te  \i/2 


X  — - — r  X  1.78  X  10^  s-\ 

27  X  10-20  m2 


Since,  in  general,  Cj  agn,  the  contribution  of  ionization  collisions  to  Ven  is  small. 


2.3.4  Electron-ion  collisions 

The  electron-ion  collision  frequency  is  expressed  as[15,  16] 


In  A  /  Ug 
12TV^/^elmy^Tp  UOi^m-^ 


1  eV\3/2 


)^(^) 


X  InA  X  2.9  X  10®  s“b 


(2.34) 


,  ,  ^  1 ,  r /10^®m  /  Ig  1 

In  A  ~  9  H —  In  f - )  ( - )  . 

2  LV  Ug  Ml  eV/  J 


Tg  \3 


Momentum  exchange  for  ions  and  electrons  is 


Mgi  =  Evgimeneivi  -  Vg). 


(2.35) 
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2.3.5  Anomalous  electron  diffusion 

Near  the  channel  exhaust,  where  the  neutral  density  is  low,  classical  collisions,  r'e*  +  i^en,  are 
insufficient  to  explain  the  plasma  behavior  observed  there,  which  suggests  the  presence  of  an 
anomalous  collisionality.  There  is  controversy  about  which  is  the  phenomenon  responsible 
for  that  collisionality.  Some  authors[17,  11]  suggest  electron  exchanges  at  lateral  walls,  but 
experimental  data[9,  18]  seems  to  favor  Bohm  diffusion.  Anomalous  diffusion  due  to  electron 
interaction  with  lateral  walls  will  be  discussed  in  Chap.  5. 

Bohm  diffusion  is  usually  included  in  electron  momentum  equation  by  adding  to  the  fre¬ 
quency  for  classical  collisions  the  Bohm  diffusion  frequency 


with 


eB 

Ug  =  —  ~ 

rrif 


^Bohm 

B 


X  1.758  X  10^  s 


9  ^-1 


0.01  Tesla 

the  electron  gyro  frequency,  and  the  Bohm  parameter,  with  a  classical  empirical  value  of 
aB  ~  1/16.  The  Bohm  diffusion  frequency  is  computed  from 

B 


^Bohm 


0.01  Tesla 


X  Was  X  1.1  X  10  s 


8  ^-1 


2.4  Equations  for  electrons 

2.4.1  Electron  current 

From  Eq.  (2.5),  the  electron  velocity  field  follows 

Dv 

meUg^^  “  -VueTe  -|-  eUgV  (j)  -  engVg  AB  +  M'^  +  (2.36) 

where 

D/Dt  =  d/dt  +  Ve  -V 

is  the  substantial  derivative. 

Mg  ~  Mg  ~  -mgng{vgVg  -  VgiVi)  (2.37) 

is  the  frictional  force,  and  the  collision  frequency 

—  ^en  T  T^gi  T  l^Bohm-  (2.38) 

includes  the  different  collisional  processes;  in  Eq.  (2.37)  we  assumed  r'j  Vg  and  J^eVg. 

The  conditions  for  a  quasi-closed,  diffusive  motion  of  the  electron  fluid  are 

D / Dt  Ug  ujg.  (2.39) 

The  condition  on  the  left,  allows  us  to  neglect  inertia  effects  in  Eq.  (2.36),  which  becomes 
algebraic  for  the  velocity  field, 

0  ~  -VugTg  -F  erieV^  -  eUgVg  AB  +  mgng{-VgVg  -F  VgiVi). 


(2.40) 
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This  simplifies  the  problem  greatly,  since  writing 


with 


0  ~  T’  —  UeVe  A  eB  —  meVeneVe, 


F  =  eneS/(p  -  ^ngTe  +  VeimeUeVi 


(2.41) 


the  total  driving  force,  we  can  solve  for  the  velocity  field.  Using  the  right  condition  in  Eq.(2.39), 
the  velocity  components  satisfy 


0  —  T 

-  F^  =  meUeiuJeVee  “  T^eVxe)  -  UJemeUeVge, 
F^ 


(2.42) 

(2.43) 

(2.44) 


where  Ug  is  the  electron  gyro-frequency.  These  equations  show  that  the  dominant  electron 
motion  is  azimuthal  and  is  due  the  F  A  B  drift,  whereas  the  axial  motion  is  due  to  the  drift 
generated  by  the  magnetic  field  and  the  azimuthal  frictional  force.  Solving  Eqs.  (2.42)- (2.43) 
for  the  axial  component  one  has 


f^e'^xe 


with 


— 


(2.45) 


(2.46) 


the  effective  collision  frequency  for  the  electron  axial  motion;  notice  that  VeimeUeVi  can  be 
neglected  in  F^. 

More  adequate  than  Vd  to  evaluate  the  relevance  of  electron  diffusion  in  the  whole  plasma 
response  is  going  to  be  the  frequency 


rUp 


rrij 


i^p 


(2.47) 


where 


eB 


B 


CJr  = 


X  3.585  X  10*^  s 


6  ^-1 


^yrrUml  0.01  Tesla 
is  the  electron  hybrid  gyrofrequency.  Estimates  of  i>d  can  be  obtained  from 

B  \2/10Vi' 


( - X  1-29  X  10®  s-^ 

VO.Ol  Tesla/  \  Vg  ) 


(2.48) 


2.4.2  Internal  energy  losses 

Taking  into  account  that  i>i  Vg  and  rUg  rrii,  the  term  for  internal  energy  losses,  Eq.  (2.8), 
can  be  reduced  to 


Q'e  —  ~'^e  ■  +  {Qi  +  Qgxc)- 


(2.49) 
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The  first  term  on  the  right  is  the  Joule  heating,  which  using  Eqs.(2.43)  and  (2.46),  satisfies 

-Ve  ■  Mg  ~  meUeil'eVl  -  l^eiVe  '  Vi)  ~  menel^eVge  -  menel^dvle'  (2.50) 

Ionization  losses  are 

Qi  —  SiEij 

and  Qexc  accounts  for  excitation  energies  that  are  lost  by  radiation.  It  is  common  to  group 
ionization  and  excitation  losses  in  a  single  term  proportional  to  Qf 

Qi  T  Qexc  —  (^iQi  —  ^if^e^^i^it 


with  Q!j  given  by 

^  1 

ai{TQ  ~  2  +  -  exp— , 

according  to  Dugan  et  o/.[19].  [However,  in  most  of  this  work  we  will  just  take  =  const  =  2.5.] 


2.4.3  Heat  conduction 

To  close  the  equations  for  electron  dynamics,  a  transport  equation  is  needed  for  the  heat 
conduction  term  V  •  q^.  The  diffusive  behavior  of  the  electrons,  expressed  by  condition  (2.39) 
justifies  the  following  diffusive  law  for  [20] 

5 

-PeVTe  +  eq^AB  +  meVeq^  ~  0,  (2.51) 


which  can  be  expressed  as 


Qe  H - Qe 

^P. 


-KoVTg,  Ko  = 


5neTg 

2mgVg  ■ 


The  azimuthal  and  axial  components  of  q^  satisfy 


0 


~^eQxe  ^eQOe 


5pe  dTg 

0  ^  - - +UJeqee, 

zrUg  ox 

where  from  we  have  the  following  Fourier-like  law  for  the  axial  component. 


Qxe 


bngTg  dTg 

2mei>d  dx 


(2.52) 


(2.53) 


(2.54) 
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2.5  Summary  of  axial  equations 

Taking  into  account  all  the  effects  discussed  above,  the  macroscopic  equations  for  the  axial 
dynamics  of  a  weakly  divergent  jet  of  quasineutral  plasma  {ug  ~  rij)  are 

OtIq  1  0 

OtIq  1  0 

dt  A  dx^  ^  e  ,  T  w, 

+  -T-^{AnnVn)  =  -UgVi  -  S-^, 

at  A  ox 

0  15 

—  {miTlnVn)  +  —  —  {AmiTlnVl)  =  -miVnUgl^i  +  mi{v^i  -  Vn)neUc^  -  miVnvjSyj, 

ot  A  ox 

d  \  ^  ^  (  A  2\  /  \  o 

—  {mingV^i)  +  —  —  {AmingV^i)  =  -(^^g—  +  rUiVnUgUi  -  mAv^i  -  VnjUgUg^  +  rriiV^iS^, 

n  9<pd 

0  =  eug- - -^ngTg  -  UdmgUgV^g, 

ox  ox 


5/3 


)  ^ (UcneV„  +  qA  =  n.  “  Te 


dVxt 

dx 


—  ViajK 


)  +  Q'we^ 


dTg  _  2mgVd^ 
dx  hUgTg 


(2.55) 


The  ion  temperature  and  pressure  have  been  neglected,  based  on  the  consideration  that  the 
region  of  ion  production  is  thin  and  therefore  Tj  Tg.  Inside  the  channel  A{x)  is  constant 
and  drops  off  from  Eqs.  (2.55),  and  the  law  A{x)  for  the  divergent  area  of  the  plasma  plume  is 
derived  next. 

2.5.1  Plume  divergence  law 

The  external  plasma  plume  is  considered  a  cylindrically  divergent  jet.  with  the  same  angle  of 
divergence,  5,  at  the  inner  and  outer  boundaries.  Thus  the  radial  area,  A(^x),  and  thickness, 
d{x)  =  r2{x)  —  ri(a;),  satisfy 


d  5  2 

—  In  ^  =  —  In  a  =  -  tan  o, 
dx  dx  a 


(2.56) 


The  simplest  approximation [21]  would  be  to  assume  6  =  const.  Instead,  experimental  mea¬ 
surements  by  Pollard  and  Beiting[22]  indicate  that  what  remains  practically  constant  (in  the 
near-plume)  is  the  radial  velocity  of  expansion  of  the  plume  boundaries,  which  we  will  call  Cp. 
This  suggests  a  law 


tan  5  =  — . 

'^xi 


(2.57) 
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The  value  for  Cp  will  be  discussed  later  for  each  particular  model. 

Of  course,  this  1-D  model  of  the  plume  will  be  more  correct  the  smaller  is  the  divergent 
angle,  that  is  the  largest  is  the  relative  ion  velocity  at  the  thruster  exhaust.  Also,  we  expect  the 
model  to  be  correct  only  in  the  near-plume;  far-plume  profiles  are  computed  with  this  model 
just  to  estimate  the  thrust  and  the  thrust  efficiency. 

2.6  The  anode  sheath 

The  Debye  sheath  attached  to  the  anode  (region  AB  in  Fig.  2.1)  completes  the  1-D  model  of 
the  channel  and  defines  two  boundary  conditions  at  point  B  for  the  quasineutral  model.  The 
sheath  is  needed  to  assure  the  continuity  of  the  electric  current  from  the  channel  to  the  anode. 
A  space-charge  field  adjusts  the  potential  jump  in  the  sheath,  =  4^b  —  <Pa-,  to  a  value  such 
that  the  flow  of  electrons  reaching  the  anode  is  equal  to  the  diffusive  flow  coming  from  the 
quasineutral  channel. 

In  the  distinguished  limit  Xd/L  -A’  0,  the  problem  in  the  sheath  is  quasisteady  and  colli¬ 
sionless.  Conservation  equations  across  the  sheath  are: 

UiV^i  =  const  =  {niV^i)B, 

\mivL  ^e(t)  =  const  =  -F  ecp^h, 

^e^xe  const  {j^e^xe) At  ^2  58) 

Te  In  Ug  —  e(f)  =  const  =  TgB  In  ngB  —  ^4^8, 


5  5 

neVxe{-^Tg  -  ecf)  +  =  const  =  {-UgV^gTg  +  q^g)A. 

Assuming  a  quasi-Maxwellian  distribution  with  temperature  Te^,  the  diffusive  fluxes  of  particles 
and  energy  deposited  by  electrons  at  the  anode  are 


/  \  ^4^sh 

[neVxejA  = - OXp  — - , 

4  leB 

{^xe  T  —  ‘^'I'esiEe'^xe) At 


(2.59) 


with  CeB  the  electron  thermal  velocity  at  B. 

Then,  the  potential  jump  across  the  sheath,  and  the  heat  conduction  at  the  sheath 
boundary,  qxeBt  satisfy 


,  CgB 

- =  fn  — - j-, 

TeB 

(2.60) 

QxeB  ^4^sh  1 

BVxeBTeB  T^b  2 

(2,61) 
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The  deposition  of  electrons  and  ions  at  the  anode  produces  anode  heating,  which  should  be 
considered  in  the  evaluation  of  thruster  performances.  The  heat  flow  deposited  by  the  electrons 
at  the  anode  is 


QeA  —  2(T'g?T.e  Vxe\)B2^c- 


Back-flowing  ions  are  accelerated  by  the  sheath  and  impinge  the  anode  with  a  large  kinetic 
energy.  In  the  anode  they  are  recombined  and  then  are  re-emitted  to  the  channel  with  a  much 
lower  kinetic  energy,  of  the  order  of  the  anode  temperature.  The  heat  flow  deposited  by  the 
ions  is 


QiA  ^  -mivliAneB\vxiB\A  =  {-mivliB  +  e(t>sh)neB\VxiB\Ac. 

Thus,  the  total  heat  flow  to  the  anode  is 


Qa  —  QiA  T  QeA  —  TgB'IT'eB^c 


‘^I'^xcbI  +  (7  + 
0 


^eB 

M'^xeB 


)  I  '^xiB  I  • 


(2.62) 


(2.63) 


The  internal  structure  of  the  sheath  is  obtained  from  the  integration  of  Poisson  equation 
with  the  aid  of  the  above  conservation  equations.  Boundary  conditions  for  Poisson  equation 
are: 

i)  in  the  sheath  scale,  the  electric  field  is  zero  at  the  presheath/sheath  transition. 


d{ecf>/TQ 


d{x/XD) 

ii)  <I)A  =  0. 

A  first  integral  of  Poisson  equation  is 

Co  /#^  2 


0. 


where. 


U  {e<j))  =  UeTe  +  miTlivli 


(2.64) 


(2.65) 

(2.66) 


is  known  as  the  Sagdeev’s  potential  and  condition  (2.64)  has  been  applied.  At  point  B,  the 
quasineutrality  condition  yields 


U{e(j))'B  =  neB  -  riiB  =  0. 


Therefore,  for  a  valid  solution  to  exist  around  point  B,  Eq.  (2.65)  requires 


U''  = 


f^eB 


n-iB 


TeB  rriiU^ 


>0, 


iB 


that  is 


2  ^ 

"^xiB  ^ 


rrij 


(2.67) 


which  is  known  at  Bohm  condition.  Conditions  (2.60),  (2.61),  and  (2.67)  are  boundary  condi¬ 
tions  for  the  quasineutral  plasma  at  point  B. 

Finally,  we  point  out  that  a  1-D  model  of  the  sheath  and  the  region  around  it  is  an  idealized 
situation,  since,  in  general,  the  gas  is  ejected  from  a  small  orifice  into  the  channel,  yielding  a 
two-dimensional  potential  distribution. 
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Figure  2.1:  Sketches  of  (a)  the  Hall  thruster  and  (b)  the  ID  macroscopic  model.  B{x)  is  the 
magnetic  field  profile;  Vd  is  the  discharge  voltage;  Id  the  discharge  current;  m  the  fiow  of  neutral 
gas;  Fq  =  UaVaA,  {a  =  i,e,n, ...)  are  particle  flows  of  the  different  species.  Surface  P  is  the 
cathode  (beam  neutralizer).  For  usual  operating  conditions  the  following  plasma  regions  will  be 
identified:  anode  sheath  (AB),  anode  presheath  (BC),  diffusion  region  (CD),  ionization  layer 
(DH),  internal  acceleration  region  (HE),  near-plume  (EP),  and  far-plume  (Poo).  Point  B  is  a 
singular  sonic  point  for  the  reverse  ion  fiow  and  point  S  is  a  regular  sonic  point  of  the  forward 
ion  fiow. 
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Figure  2.2:  Individual  {ai)  and  collective  (dj)  cross-section  for  single-ionization  of  xenon  by 
electron  impact;  E  is  the  energy  of  an  individual  electron  and  T^  is  the  temperature  of  the 
Maxwellian  electron  distribution.  Asterisks  and  circles  are  experimental  values  taken  from 
Refs.  [23]  and  [24],  respectively.  The  dashed  line  for  dj  corresponds  to  the  integration  of  the 
(average)  experimental  points  of  ai{E).  Solid  lines  represent  the  analytical  fittings  given  by 
Eqs.(2.27)  and  (2.28)  for  =  5  x  IQ-^W. 


Chapter  3 

BASIC  STEADY-STATE  SOLUTION 


3.1  Governing  equations 


We  make  here  the  following  simplifications  on  the  axial  model: 

(i)  steady-state  solution:  d/dt  =  Q 

(ii)  no  wall  losses:  =  0 

(iii)  no  ion-electron  collisions:  t'e*  =  0 

(iv)  no  charge-exchange  collisions:  Vcx 

(v)  no  heat  conduction  Qg  =  0 

(vi)  The  axial  profile  of  the  magnetic  field  is  taken  Gaussian, 

B  =  Bcr,  B{x)  =  -B^exp  -  , 

with  Xrm  and  constants,  and  x^  close  to  xe- 
Then,  Eqs.(2.55)  simplify  into 

^{AUeV^i)  =  -^{ArieV^e)  =  =  AUel^i, 

Vn  =  const. 


(3.1) 


dVxi  d(f)  ,  , 

^  ~^~dx  ~ 


(3.2) 


0 


tIqTq  z^^rngUgii^g, 

ax 


^  ^  +  ^Tg). 

These  equations  are  completed  with  Eq.  (2.57)  for  the  variation  of  the  radial  area,  a  =  Ind  = 
In  A  -h  const,  in  the  plume. 

Using  the  particle  axial  flows  of  the  different  species  as  Fq  =  AuaVxa,  the  mass  flow  and 
the  discharge  current  are  written  as 


771  —  TJljVjn^  —  U  -h  F^, 

Id  =  eFrf,  Frf  =  Fj  —  Fg. 
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Then,  continuity  equations  show  that  these  two  magnitudes  are  constant  along  the  plasma  jet 
(except  for  dropping  abruptly  to  zero  at  the  cathode).  Notice  that  m  includes  only  the  gas 
flow  delivered  at  the  anode.  Another  useful  equation,  derived  from  Eqs.(3.2),  is  the  ambipolar 
momentum  equation 


d 

dx 


{AmiUevli)  +  A-^{neTe) 


ArieiVirriiVn  -  VdmeVxe)- 


(3.3) 


Finally,  important  ratios  for  a  discussion  of  the  plasma  response  are  (i)  the  ratio  between 
the  discharge  current  and  the  emitted  gas  flow. 


id 


rrijld  ^ 
em  “  r„’ 


and  (ii)  the  local  ion  current  fraction. 


r]i{x)  =  Vi{x)/V^ 


(3.4) 


(3.5) 


The  propellant  utilization  is  defined  as  For  id  >  I  (the  usual  case  in  practice)  can 

be  near  unity[17];  in  this  case,  coincides  with  the  electrical  efficiency  However,  for 

id  <  1,  the  maximum  propellant  utilization  can  be  just  id- 


3.1.1  Singular/sonic  points 

We  discuss  now  the  properties  of  the  mathematical  model.  Solving  Eqs.(3.2)  for  the  spatial 
derivatives  we  obtain  a  matrix  equation  of  the  form 

dY 

{1-M^)—  =  F{Y),  (3.6) 

where  Y  =  {ne,Vxi,Te, ...)  groups  7  plasma  variables,  E’  is  a  regular  function  (different  from 
the  one  in  the  zero  conduction  model),  and  M  is  the  isentropic  Mach  number. 


M 


(3,7) 


For  instance,  the  equations  for  rig  and  are 


with 


dUf. 

dx 

dVxi 

dx 


n. 


G 

e  p  1 


Vi 


G  da\ 
P  ^ 


(3.8) 


r2 


P  =  Tg--m,nT  =  Tg(l-M2), 


Q  3  Q  dOi 


(3.9) 
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Therefore,  the  spatial  problem  may  present  singular  points  where  the  Mach  number  is  equal 
to  -1  or  +1.  The  definition  of  an  ion  Mach  number  is  natural  to  the  macroscopic  approach 
adopted  to  describe  plasma  dynamics.  It  characterizes  the  propagation  of  pressure  perturba¬ 
tions,  and  clarifies,  therefore,  the  interpretation  of  the  steady  and  unsteady  plasma  responses. 
Sonic  points  are  regular  or  singular  depending  on  G  being  equal  or  not  to  zero  at  them.  In  a 
singular  sonic  point  the  derivatives  of  all  plasma  variables  become  infinite,  so  it  is  a  turning 
point  of  spatial  profiles.  These  points  are  admissible  only  at  the  boundaries  of  the  quasineutral 
channel.  On  the  contrary,  in  a  regular  sonic  point  the  ion  flow  changes  smoothly  from  subsonic 
to  supersonic,  like  in  the  classical  problem  of  a  convergent-divergent  nozzle.  The  expansion  of 
the  ion  beam  from  the  channel  into  the  rarefied  ionosphere  implies  that  the  ion  flow  will  be 
sonic  or  supersonic  at  the  channel  exit  (point  E).  Therefore,  a  forward  sonic  point  -point  S 
in  Fig.  2.1(b)-  with  Ms  =  1  is  always  expected  in  the  channel.  For  a  choked  exit,  point  S  is 
singular  and  coincides  with  point  E.  For  a  supersonic  exit,  point  S  is  regular  and  its  position 
within  the  channel  must  be  determined  as  part  of  the  solution. 

3.1.2  Boundary  conditions 

The  integration  of  Eqs.  (3.2)  between  points  B  and  P  require  seven  boundary  conditions,  which 
are  distributed  between  points  B,  S,  E,  and  P“.  Taking  into  consideration  the  singular  points, 
the  sheath  conditions  discussed  in  Sec. 2. 6  and  the  parameters  that  are  controlled  experimentally, 
a  natural  set  of  boundary  conditions  is  the  following 

i)-ii)  The  externally  supplied  mass  flow  rate  m  and  the  neutral  velocity  VnA  are  known. 

iii)  The  electron  temperature  at  the  cathode,  T^p^  (of  the  order  of  the  electron  emission 
temperature)  is  known. 

iv)  The  discharge  voltage  =  (f)A  —  (t>p  known. 

v)  The  quasineutral  solution  in  the  channel  requires  that  —  <  1  [were  —Mp  >  1,  there 

would  be  a  point  within  the  channel  with  M  =  —1  and  G  =  0,  which  is  not  possible,  as 
the  solution  will  make  evident.]  On  the  other  hand,  the  Bohm  condition  (2.67)  at  the  sheath 
entrance  states  that  —Mp  >  -y/S/S.  Here,  we  will  assume  that  the  sheath  transition  occurs  at 
the  singularity  of  the  quasineutral  solution,  that  is 

Mp  =  -1.  (3.10) 

vi)  The  voltage  difference  at  the  sheath  satisfies  Eq.  (2.60).  (Placing  the  arbitrary  origin  of 
(j)  at  point  A,  one  has  (j)p  =  (psh  and  (l)p  =  —Vd).  Using  Eq.(3.10),  Eq.(2.60)  becomes 

TeB  ^  VPeB  V  lOTTme/ 

vii)  For  the  supersonic  exit  regime^  sonic  point  S  is  regular: 

Gs  =  0  at  Ms  =  1.  (3.12) 

Instead,  for  the  choked  exit  regime  we  set  that  the  flow  at  point  E  is  sonic: 


(3.11) 


Mp  —  1. 


(3.13) 
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For  the  velocity  of  radial  expansion  in  the  plume,  Eq.  (2.57),  we  make  the  choice 


Cp  — 


(3.14) 


Apart  from  V^,  m,  and  Tgp,  other  input  variables  are  v^a,  L,  Lep,  A^.  and  the  profile  of  the 
magnetic  field,  Eq.(3.1).  The  main  output  variables  of  the  model  are  the  discharge  current  I^, 
the  relative  ion  reverse  flow  rjis,  the  position  xs  of  the  sonic  point  S  (for  a  supersonic  exit),  the 
plume  divergence  profile,  the  thrust,  and  the  thrust  efficiency. 

The  continuation  of  the  solution  from  point  P  to  oo  is  straightforward.  It  just  requires  us 
to  modify  the  electron  current  at  point  P+,  setting  P+p  =  Pjp,  in  order  that  the  net  electric 
current  in  the  plume  far-held  be  zero. 


3.2  Types  of  solutions 

It  is  convenient  to  use  dimensionless  variables  and  parameters  to  integrate  the  equations  and  to 
characterize  the  solutions.  Three  magnitudes  are  enough  to  nondimensionalize  Eqs.(3.2).  The 
energy  T^,  =  Ei,  the  cross-section  a*  =  ^  and  a  typical  particle  how  P*  are  the 

best  choice.  Then,  magnitudes  related  to  these  three  are 

—  '\Y  T*  5  —  P*/Acr)j(!,  U  —  —  v^jU. 

When  necessary,  dimensionless  magnitudes  will  be  represented  by  an  over-tilde  over  the  respec¬ 
tive  dimensional  ones.  We  omit  to  write  down  the  dimensionless  form  of  plasma  equations, 
which  is  immediate. 

The  dimensionless  form  of  Eqs.  (3.2)  is  integrated  as  an  initial- value  problem  with  a  Runge- 
Kutta  algorithm.  We  focus  the  discussion  here  in  the  region  between  points  B  and  P,  the 
continuation  to  the  far-plume  being  straightforward.  The  use  of  the  auxiliary  variable  dehned 
by 

dx/d^  =  P, 

instead  of  x,  avoids  unnecessary  numerical  problems  in  reaching  or  departing  from  singular 
sonic  points,  like  point  B.  However,  this  does  not  solve  the  crossing  of  regular  sonic  point  S 
in  the  supersonic-exit  regime.  A  Taylor  expansion  of  Eqs.  (3.2)  is  needed  to  determine  the 
derivatives  of  the  plasma  variables  there;  the  details  of  that  expansion  are  omitted  here. 

Due  to  the  character  of  point  S,  the  less  time-consuming  procedure  is  clearly  to  integrate 
from  the  neighborhood  of  point  S  towards  points  B  and  P  separately.  At  point  B,  the  integration 
stops  where  condition  Mb  =  —1  is  reached;  at  point  P,  it  ends  where  T  =  Tep-  This  initial- 
value  procedure  requires  setting  all  variables  at  point  S.  For  the  choked-exit  regime,  this  implies 
taking  output  parameters  rjis,  id,  and  T^s  as  initial  parameters,  while  input  parameters  V^,  xa, 
and  xp  are  outputs  of  the  integration.  For  the  supersonic-exit  regime,  output  parameters 
r/js,  and  xs  are  used  as  initial  parameters,  and  TgS  comes  from  condition  Gs  =  0;  as  before, 
input  parameters  xa,  and  xp  are  outputs  of  the  integration.  A  subsequent  iteration  is 
needed  to  adjust  these  input  parameters  to  their  desired  values. 

In  addition,  not  all  sets  of  ’initial  parameters’  lead  to  valid  solutions.  A  first  restriction 
to  these  parameters  comes  from  the  condition  that  the  plasma  flow  be  accelerating  at  point 
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S:  dvxi/dx\s  >  0.  For  a  supersonic-exit  solution  this  means  choosing  the  appropriate  set  of 
derivatives  from  the  two  solutions  given  by  the  Taylor  expansion.  For  a  choked-exit  solution  it 
reduces  the  region  of  valid  parameters  to  those  satisfying 


>  0  and  Gg  <  0; 


(3.15) 


observe  in  Eq.  (3.9)  that  the  difference  between  Gg  and  Gg  is  due  to  the  term  containing 
da/dx  oc  (tan5)/d. 

A  second  restriction  to  ’initial  parameters’  at  point  S  comes  from  the  fact  that  there  are 
solutions  departing  from  S“  towards  B,  which  never  reach  a  negative  ion  velocity  and  finish 
instead  in  a  singular  point  with  M  =  -hi.  These  solutions  are  disregarded  since  they  have  no 
physical  meaning  for  the  plasma  in  a  conventional  Hall  thruster.  The  limit  condition  Fj^  ~  0 
is  going  to  be  very  important  since  it  is  going  to  indicate  a  real  physical  condition  for  the 
disappearance  of  steady  solutions. 

Figure  3.1  illustrates,  in  the  parametric  plane  {Tes,Vis)  and  for  two  values  of  id,  Eq.(3.4), 
the  different  solutions  that  are  obtained  from  the  integration  between  points  S  and  B.  Four 
parametric  regions  are  found:  region  1,  with  G~^  <  0,  and  region  2,  with  Gg  >  0  and  rjiB  >  0, 
are  forbidden  regions  for  the  reasons  exposed  above;  region  3,  with  G~^  >  0  and  rjis  <  0, 
corresponds  to  the  choked-exit  regime;  and  (line)  region  4,  with  Gs  =  0  and  rjis  <  0,  corresponds 
to  the  supersonic-exit  regime.  For  the  choked-exit  regime,  the  continuity  of  the  solutions  into 
the  plume  requires  that  G^  <  0,  Eq.(3.15).  Figure  3.1  depicts  the  limit  line  G^  =  0  for  two 
values  of  (tan5)/d. 

For  id  =  1.25,  region  4  of  supersonic-exit  solutions  consists  of  high  and  low  ionization 
branches.  The  low-ionization  branch  yields  propellant  utilizations  well  below  50%  and  corre¬ 
sponds,  then,  to  inefficient  thruster  operation.  Since,  in  addition,  it  does  not  lead  to  different 
plasma  dynamics  than  the  high-ionization  branch,  it  will  be  disregarded  hereafter  in  this  pa¬ 
per.  The  high-ionization  branch  is  actually  a  near-total  ionization  branch  and  corresponds  to 
the  near-complete  ionization  mode  referred  by  Morozov  et  o/.[17].  Notice  that,  since  rjis  —  1, 
point  S  must  be  placed  near  the  downstream  end  of  the  ionization  layer  (point  H  in  Fig.  2.1). 
As  id  increases  from  1.25,  the  line  rjis  —  0  (separating  regions  3  and  4  from  region  2)  moves 
left,  towards  lower  T^s,  and  region  2  disappears  eventually  (at  id  ~  1.65  for  the  case  of  Fig. 
3.1).  For  the  same  reason,  as  id  decreases  from  1.25  the  line  rjis  —  0  moves  right  and  restricts, 
in  particular,  the  domain  of  supersonic-exit  solutions  with  large  ionization.  For  <  1,  two 
features  are  to  be  pointed  out:  first,  the  high-ionization  branch  of  the  supersonic-exit  regime, 
instead  of  approaching  the  asymptotic  line  Fj^  =  F,^  (i.e.,  rjis  =  id),  disappears;  second,  the 
line  rjiB  —  0  forbids  choked-exit  solutions  with  large  ionization  and  moderate  temperatures. 

The  presence  and  characteristics  of  the  two  supersonic-exit  branches  can  be  explained  from 
the  analytical  expression  of  Gs  =  0.  From  Eq.(3.9),  that  condition  yields 


3F 


iS_ 

m 


{I'm 


p  (2F(j  —  Fj5')Te5'  -h  {2/5)TisaiEi  ^ 

“  -  Os)"  ~  I-* 


(3.16) 


with  Vd  =  and  pt  =  nnARi{Tes),  and  this  last  magnitude  is  of  the  order  of  the 

maximum  ionization  frequency  in  the  channel.  Then,  for  plasmas  satisfying 


(3.17) 
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it  is  readily  seen  that  Eq.(3.16)  is  verified  by  two  types  of  solutions:  the  low-ionization  branch 
with  <C  F^;  and  the  near-total  ionization  branch  with  Fj^  ~  F^,  which  exists  only  for 
Trf  >  F^. 

In  conclusion,  we  have  shown  that  a  regular  sonic  point  (with  no  discontinuity  in  plasma 
derivatives)  may  exist  inside  a  channel  of  constant  area  in  a  natural  way  by  compensation  of 
diffusion  and  ionization  effects.  This  implies  that  the  sonic  point  is  within  or  close  to  the 
ionization  layer.  The  position  of  the  sonic  point  is  not  set  a  priori^  it  depends  on  the  control 
parameters  and  is  part  of  the  solution  of  the  problem.  Regular  sonic  transitions  were  also 
discussed  by  Fruchtman  and  Fisch[10].  In  addition,  they  suggested  forcing  a  choked  internal 
point  by  introducing  a  current  discontinuity;  thus,  their  solution  matches  two  choked  flows 
(with  Gg  >  0,  G~^  <  0)  at  the  prescribed  discontinuity  point. 


3.3  Structure  of  the  plasma  flow 

The  preceding  section  has  shown  that  supersonic-exit  solutions  with  moderately  above  1 
seems  to  be  the  most  efficient  operating  regime  of  a  Hall  thruster.  The  results  we  show  next  on 
the  plasma  structure  and  the  thruster  output  parameters  will  confirm  this.  These  results  are 
presented  in  dimensional  form  to  facilitate  the  comparison  with  real  experiments. 

Figures  3.2(a)-(i)  show  spatial  profiles  of  the  main  plasma  variables  for  a  supersonic-exit 
solution  and  geometrical  and  operational  parameters  typical  of  a  SPT-100  thruster[25];  all 
parameters  are  listed  in  the  figure  caption.  We  observe  that  the  structure  of  the  plasma  response 
can  be  divided  into  the  six  regions  sketched  in  Fig.  2.1(b).  The  first  significant  feature  is  that 
the  ionization  region  is  placed  in  the  middle  of  the  channel,  more  or  less  between  points  D  and 
H,  which  can  be  defined  as  the  places  where  v^i  changes  sign  and  Tg  is  maximum,  respectively; 
this  yields  xd  ~  13.6  mm  and  xh  —  17.4  mm,  with  T^h  ~  91.3  eV.  Notice  in  Fig.  3.2(b)  the 
peaked  shape  of  i'i{x).  The  point  of  maximum  electron  pressure,  p^  =  rigTe,  is  at  a;  ~  14.1mm, 
very  near  to  the  proposed  entrance  of  the  ionization  layer.  The  sonic  point  S  is  downstream 
of  point  H  [xs  —  19  mm)  and  the  temperature  there  {Tes  —  89.6  eV)  is  close  to  the  peak 
temperature;  notice  that  point  S  cannot  be  far  from  point  H  since  the  ionization  terms  are 
needed  to  define  the  crossing  condition  Gs  =  0,  Eq.  (3.16). 

Upstream  of  the  ionization  layer  is  the  large  dijfusion  region  (region  CD),  characterized  by 
a  reverse,  low- velocity  ion  flow,  a  low  electron  temperature  (~  2eV),  an  insignificant  electric 
field  {^c  —  ^  I  and  a  quasi-linear  pressure  gradient  as  the  force  driving  the  electron 

diffusion.  The  relative  value  of  the  reverse  ion  flow  is  \r]iB\  —  0.24  (of  the  order  of  the  one 
measured  in  Ref.  [1]).  Tied  to  the  anode,  there  is  the  anode  sheath  (region  AB)  characterized 
by  a  potential  jump  —  3.8  V  and  a  negligible  thickness  in  the  channel  scale  {xb  —  xa  =  0 
for  Xd/L  0.).  The  anode  presheath  (region  BC)  is  a  transition  region  between  the  sheath 
and  the  diffusion  region,  where  the  reverse  ion  flow  is  accelerated  to  meet  the  Bohm  condition 
at  the  sheath  entrance;  it  is  rather  thin  [xc  —  xb  ^  1mm)  and  observable  in  Fig.  3.2(f)  only. 

The  sharp  transition  from  the  diffusion  region  to  the  ionization  region  is  related  to  the 
exponential  dependence  of  the  ionization  cross-section  on  the  electron  temperature  Tg,  as  we  will 
further  comment  on  below.  The  large  value  of  Tg  in  the  ionization  region  -  and  its  low  value  in 
the  diffusion  region-  are  consequences  of  excluding  the  heat  conduction  and  the  losses  to  lateral 
walls.  The  electric  potential  remains  practically  constant  until  well  inside  the  ionization  layer. 
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then  it  starts  decreasing;  the  total  potential  variation  across  the  ionization  layer  is  (pD  —  (1>h  — 
51.4  V. 

Downstream  of  the  ionization  layer  there  is  the  internal  acceleration  region  (region  HE), 
where  the  ion  flow  is  accelerated  by  the  electric  held.  The  potential  decrement  in  this  region  is 
(pH  —  (pE  —  91.0  eV  and  the  temperature  at  the  channel  exhaust  is  still  large  {T^e  —  65.3eV)  so 
the  effects  of  the  electron  pressure  are  significant  there  still.  The  plume  is  the  extension  of  the 
internal  acceleration  region.  The  form  of  the  plume  profile  shown  in  Fig.  3.2(b)  comes  from 
the  solution  for  A{x)  in  Eqs.(3.2).  The  semi-angle  of  divergence  is  maximum  at  the  channel 
exit,  Se  —  32.3°,  but  it  does  not  change  much  since  most  of  the  increase  of  the  ion  velocity 
took  place  inside  the  channel;  the  typical  range  of  6  in  experiments[26,  27]  is  30°  —  40°.  In 
the  near-plume  (region  EP)  the  radial  area  increases  from  Ae  =  A^.  —  45  cm^  to  Ap  ~  76 
cm^.  The  potential  variation  in  the  near-plume  is  (pE  —  (pp  —  150  V,  about  53%  of  the  total 
discharge  voltage.  The  minimum  of  the  electric  potential  at  the  cathode  (point  P)  provides  the 
electric  held  required  to  move  the  electrons  into  the  beam.  Indeed,  the  potential  variation  in 
the  far-plume,  (p^  —  pp  34.8  V,  can  be  considered  as  a  part  of  the  energy  loss  associated  to 
the  cathode. 

There  are  two  Tree’  parameters  in  the  solution  of  Fig.3.2:  Lpp  and  ap-  The  first  parameter 
seems  to  affect  weakly  the  plasma  response  except  when  it  becomes  too  small.  On  the  contrary, 
the  plasma  response  is  very  dependent  on  the  value  of  the  Bohm  parameter.  In  fact,  the 
best  fit  with  experimental  results  is  obtained  with  Wap  ~  0.2  instead  of  the  classical  value. 
Wap  =  1.  For  the  chosen  value  of  anomalous  diffusion  becomes  dominant  only  downstream 
the  ionization  region,  as  the  quick  change  in  the  derivative  of  Vd  around  point  D,  in  Fig.3.2(i), 
illustrates. 

The  plasma  structure  for  a  choked-exit  solution  is  very  similar  to  the  one  presented  here 
except  that  (i)  ionization  might  not  be  near-total,  (ii)  the  region  of  supersonic  acceleration  is 
totally  external,  and  (iii)  the  plume  divergence  is  large.  Indeed  for  practical  values  of  h/L,  the 
condition  Gg  =  0,  depicted  in  Fig.3.1,  leads  to  divergence  angles  above  45°.  Clearly,  the  1-D 
plume  model  is  hard  to  apply  there.  This  consideration  and  the  fact  that  the  experimental  data 
fits  in  general  well  with  supersonic-exit  solutions  makes  advisable  to  continue  the  investigation 
of  choked-exit  solutions. 


3.4  Thruster  performances 

Looking  again  at  the  case  of  Fig.3.2,  we  show  now  that  the  main  output  parameters  for  thruster 
performances  agree  rather  well  with  experimental  values[25].  The  discharge  current  is  G  —  4.36 
A  (versus  4.5  A  experimentally)  The  thrust,  computed  at  the  channel  exhaust  as 

F  =  mi{TiV^i)E  +  {pe  -  P<x)A,  (3.18) 

is  F  ~  92.8mN  (versus  82  mN  experimentally).  The  thrust  efficiency,  defined  in  the  usual  way 
as 


r]  =  FV2m/dFd,  (3.19) 

is  r/  ~  63%  (versus  52%  experimentally,  excluding  cathode  flow).  Energy  losses  neglected  in 
this  model  will  reduce  somehow  these  two  last  magnitudes.  It  is  interesting  to  measure  how 
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different  phenomena  affect  the  thrust  efficiency.  This  can  be  factorized  as 

V  =  vlVcVeVp,  (3.20) 

with  r]u  =  rjioo  the  propellant  utilization,  Eq.(3.5),  rjc  =  1/id  the  current  efficiency,  rje  = 
mivL^/2eVoo  the  energy  efficiency,  and  rjp  =  rjpcVpd  the  plume  efficiency,  this  last  one  split 
into  the  cathode  efficiency:  rjpc  =  V^jVd-,  and  the  divergence  efficiency:  rjpd  =  F^/F^,  with 
Fqo  =  mi{riVxi)oo-  For  the  case  of  Fig.  3.2,  one  has,  first,  that  the  propellant  utilization  and 
the  energy  conversion  are  practically  total:  ~  100%,  and  rje  ~  99.3%.  Efficiency  losses  are 

distributed  among  the  other  three  contributions:  rjc  —  88.7%,  rjpc  —  88.1%,  and  rjpd  ~  81.1%. 

The  present  model  does  not  consider  the  gas  spent  in  the  cathode  for  electron  emission, 
which  amounts  to  a  5-10%  of  the  total  gas  flow,  typically [26],  and  most  of  it  (above  90%)  is 
ejected  into  the  plume  without  ionizing.  On  the  one  hand,  this  gas  flow  contributes  as  a  loss 
to  the  thruster  efficiency;  on  the  other  hand  it  is  responsible  for  charge-exchange  and  backflow 
effects  in  the  plume. 


3.5  Asymptotic  analysis 


We  derive  here  approximate  solutions  for  the  different  plasma  regions  that  explain  most  of  the 
features  observed  in  the  exact  solution  studied  in  the  preceding  section. 

Since  ionization  is  limited  to  a  thin  layer,  the  zero-ionization  form  of  the  plasma  equations 
will  be  applicable  to  most  of  the  plasma  flow.  Making  — >■  0  in  Eqs.  (3.2),  one  has  conservation 
equations  for  particles,  ion  energy  and  electron  enthalpy, 

Tq  =  const  {a  =  i,e,n),  (3-21) 

mj'y^j/2  -\-  ecj)  =  const,  (3.22) 

NT^j^  —  e(f)  =  const.  (3.23) 


To  avoid  integral  terms  with  the  electron  diffusion  frequency  and  simplify  thus  the  analytical 
expressions,  we  will  use,  instead  of  a  continuous  function,  two  constant  (i.e.,  average)  values: 

=  ^dOi  for  the  whole  back-flow  region  AD,  and  Vd  =  for  the  whole  acceleration  region 
DP. 

We  proceed  to  integrate  Eqs. (3.2)  from  point  B  (with  ~  =  0,  Mb  =  —  1)  towards 

the  cathode  P.  To  neglect  ionization  in  the  rear  part  of  the  channel  we  assume  that  <C  Ej. 
Then,  from  Eqs.(3.21)-(3.23)  and  (3.3)  the  ion  velocity  satisfies 


'^xi 


'^xiB 


(3.24) 


where  v^w  =  —  \/5Te/3mj  and 


8  771^  VxiB  F{5 

5  mel^dO  FeB 


is  the  characteristic  thickness  of  the  anode  presheath.  Observe  that  xc  tends  to  zero  with  Pj^. 

The  diffusion  region  corresponds  to  a;  xc-  From  Eq.  (3.24),  the  asymptotic  behavior  of 
the  ion  velocity  is 

VxihxiB  —  Xc/M  -)■  0. 
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The  plasma  temperature  and  the  voltage  also  tend  to  constant  values, 

Te  — )•  Tec  =  (4/3)TeB,  e(f)  — >■  e(l)c  =  e(f)B  +  (5/6)TeB,  (3.25) 

so  the  electric  field  tends  to  zero.  On  the  contrary,  the  electron  pressure  increases  monotonically: 

Pe  -  PeB  +  l^e/AclrUeVdoX,  (3.26) 


and  provides  the  driving  force  for  the  electron  axial  motion  in  the  diffusion  region.  Since 
Te  ~  const,  the  plasma  density  changes  in  the  same  way  as  the  pressure.  Weak  electric  fields 
and  nonzero  density  gradients  upstream  of  the  ionization  layer  are  confirmed  by  experiments[l]. 

The  transition  to  the  ionization  region  (around  point  D)  requires  the  inclusion  of  small 
ionization  effects.  The  evolution  of  Tg  there  is  governed  by  the  exponential  dependence  of  the 
ionization  rate  on  Tg,  for  Tg/Tj  <C  1.  Keeping  d^/dx  ~  0  and  using  Eq.(3.26)  for  p^,  the 
electron  temperature  in  this  transition  layer  satisfies 


dTf.  Te  +  aiEi 

— ~  - — - meVdoxdx 

Ri{Te)  "  Te 

(dj  =  2Q!j/5),  and  its  approximate  solution  is 

2  _  ^2 _ ^3 _ 

{Te  +  aiEi)Eimer'donnRi{Te)' 

with 


(3.27) 


^AD 


27;c 


ajEf 


Tec 

rrieizdoi^ic ' 


(3.28) 


Equation  (3.27)  yields  Te(a;)  implicitly,  showing  a  large  gradient  around  Te/Tj  =  0(1),  where 
the  ionization  rate  Ri{Te)  increases  sharply.  This  justifies  that  Lad,  as  defined  by  the  two  last 
equations,  is  effectively  the  length  of  region  AD.  Equation  (3.28)  shows  also  how  the  position 
of  the  ionization  layer  is  influenced  by  different  factors.  The  pressure  satisfies  Eq.  (3.26)  up  to 
the  ionization  layer  and,  setting  its  maximum  around  point  D,  one  has 


PeD  ^  \TeB/Ae\meT^doLAD- 


(3.29) 


The  plasma  density,  obtained  from  Pe/Te,  gets  its  maximum  before  point  D,  due  to  the  fast 
increase  of  Tg. 

A  complete  analytical  solution  for  the  ionization  layer  is  not  available,  but  there  are  several 
useful  expressions  relating  point  D  and  point  S  (just  downstream  of  point  H,  as  we  commented 
on  previously).  Conditions  Vnjvxis  ^  1  and  (3.17)  are  enough,  generally,  to  neglect  electron 
diffusion  in  the  bulk  of  the  ionization  layer,  and  the  right-hand  side  of  the  ambipolar  momentum 
equation  (3.3).  Then,  one  has 

mnevli  +Pe  ^PeD, 

across  that  layer.  Particularizing  at  point  S  and  using  Ms  =  1  we  obtain  two  relations. 


3 

PeS  '^PeD, 


TeS  = 


3m  jT?^ 


(3.30) 
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and  ~  F^  for  a  supersonic-exit  solution.  From  Eqs.(3.2)  and  (3.29),  the  estimated  thickness 
of  the  ionization  layer  is 

r  Y  '^xiS  Y 

Ldh  ~  J-'DS  ~  ~  - r^AD, 

vl  rriiVl 


and  Lds  ^  Lad^  as  expected.  From  the  equation  for  ion  momentum  in  Eqs.(3.2),  the  potential 
jump  across  the  ionization  layer  is 


e(j)D  -  e(j)s  -Tes  +  rrii  I  v^iVidx  =  (1  +  di)-T^s,  (3.31) 

6  Jxn  6 

with  parameter  5i  small  (~  0.10)  in  the  cases  we  have  analyzed. 

In  the  acceleration  region  and  the  plume  of  a  supersonic-exit  solution,  there  are  practically 
no  neutrals  to  ionize  and  Eqs.(3.21)-(3.23)  are  applicable  again.  For  T^p  <C  and  rjis  —  1 
one  finds 


VxiP  -  e(f)s  -  e(j)p  ~  {b/2)T^s-  (3.32) 

The  area  variation  in  the  plume  frustrates  an  exact  integration  of  Eq.(3.3).  Nevertheless,  an 
approximate  integration,  useful  to  evaluate  the  influence  of  different  parameters,  is 

3(Fd  —  F^) 

PeS  ^  ^ - rUel^dlLsP, 

with  A  an  average  value  of  the  the  radial  area  of  the  jet.  This  completes  the  analysis  on  the 
different  plasma  regions. 


3.6  Influence  of  control/design  parameters 

We  restrict  the  discussion  to  supersonic-exit  solutions.  Before  showing  exact  numerical  results, 
let  us  use  the  preceding  asymptotic  analysis  on  the  plasma  behavior,  to  obtain  an  approximate 
idea  of  the  scaling  laws  among  the  the  thruster  parameters.  Since  that  analysis  corresponds  to 
a  rather  ideal  situation  (no  losses  to  lateral  walls  and  no  heat  conduction)  some  laws  must  be 
taken  with  prudence. 

First,  adding  the  potential  jumps  in  the  different  regions  [Eqs.(3.25),  (3.31),  and  (3.32)], 
one  has. 


OYd  +  Q<Psh  —  H - ^ - -TeS]  (3.33) 

observe  that  the  potential  discharge  in  the  acceleration  region  is  about  three  times  that  of  the 
ionization  region.  Next,  keeping  only  dominant  terms,  the  maximum  plasma  temperature  turns 
to  be  proportional  to  the  discharge  voltage, 

TeH  ~  Tes  ^  {3/10)eVd. 

This  law  agrees  well  with  the  numerical  results  of  Fig.  3.2,  which  indicates  that  the  overestimate 
of  the  maximum  temperature  is  intrinsic  to  the  present  model.  Therefore,  energy  losses  at 
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lateral  walls  and  heat  diffusion  cannot  be  excluded  in  order  to  estimate  correctly  the  plasma 
temperature. 

From  Eq.(3.32)  and  for  Lep  0,  the  thrust,  Eq.(3.18),  satisfies 


(3.34) 


From  Eq.(3.30),  the  maximum  pressure  in  the  channel  verifies 


8  4F 


(3.35) 


The  mean  temperature  in  the  diffusion  region  is  T^c  —  (4/3)Tes.  From  Eq.  (3.28),  one  has 


Tec  oc  In  ^  (vlVamfud) . 

From  Eq.(3.5),  the  length  of  the  acceleration  region  verifies 

T  ~  1  [Wd 

“  (1  -  r]c)  V  m  ’ 

which  differs  from  the  law 

Use  ^  VQ^d/rUiUdUi, 

proposed  by  Erofeev  and  Zharinov  (cited  in  [5]).  Adding  Lse  to  Lad,  deduced  from  Eq.  (3.29), 
one  has 


(3.36) 

(3.37) 


4r?c  ^ 


4(1  -  r/c)z^dij  ’ 


which  is  an  implicit  equation  for  r/c. 


(3.38) 


Figures  3.3  to  3.8  show  the  effects  of  V^,  m,  5^,  L,  Lep,  and  ap  on  the  plasma  response 
and  the  thruster  performances.  These  plots  have  been  obtained  from  continuation  routines, 
departing  from  the  solution  of  Fig.  3.2.  The  main  remarks  are: 

1.  The  thrust  depends  on  the  gas  flow  and  the  discharge  voltage,  and  follows  rather  closely 
Eq.(3.34),  which  is  in  good  agreement  with  experimental  data. 

2.  Contrary  to  the  thrust, the  efficiency  rj  depend  on  (and  increase  with)  the  channel  length 
and  the  magnetic  field.  The  increase  with  agrees  well  with  experimental  data  [5]. 
However,  experiments  indicate  often  that  rj  presents  a  maximum  with  L,  instead  of  the 
monotonic  behavior  found  here;  losses  at  lateral  walls,  which  seem  to  increase  with  L, 
would  explain  the  difference. 

3.  Since  the  energy  and  ionization  efficiencies  are  practically  total  (r/„,  rje  ~  1)  in  our  model, 
the  variation  of  the  thruster  efficiency  depends  basically  on  rjc  and  rjp;  notice  also  that  rjc 
is  a  dimensionless  representation  of  the  inverse  of  the  discharge  current,  1/Id- 


CHAPTER  3.  BASIC  STEADY-STATE  SOLUTION 


28 


4.  Contrary  to  common  ideas,  the  distance  from  the  ionization  layer  to  the  thruster  exhaust 
depends  more  on  and  L  than  on  the  magnetic  field  profile  and  strength.  For  usual 
conditions,  the  ionization  layer  is  placed  in  the  outer  half  of  the  channel,  and  never  next 
to  the  anode.  This  is  consequence  of  the  role  of  the  electron  pressure  on  the  plasma 
response. 

5.  The  average  temperature  in  the  inner  part  of  the  thruster,  T^c,  depends  mainly  on  m  and 
L.  In  our  model,  possibly  because  of  the  zero  heat  conduction,  this  temperature  is  too 
low. 

6.  Except  m,  the  rest  of  the  control  parameters  affect  significantly  the  relative  reverse  flow 
of  ions,  rjiB-  Experimental  data  on  the  ion  flow  in  the  deepest  part  of  the  thruster  is 
difficult  to  obtain.  Few  experiments  have  identified  the  existence  of  a  reverse  flow  and 
we  did  not  find  anyone  that  measures  the  changes  in  the  reverse  flow  with  the  control 
parameters. 

7.  There  is  an  indirect  experimental  confirmation  that  our  predictions  on  the  reverse  ion 
flow  are  correct.  Our  model  predicts  that  stationary  solutions  fail  to  exist  when  the 
plasma  cannot  sustain  a  reverse  ion  flow;  this  corresponds  to  points  with  rjis  —  0  in  the 
Figures.  If  we  take,  for  instance,  the  case  of  Fig.  3.5,  where  the  magnetic  field  strength  is 
varied,  we  observe  that  there  is  a  maximum  value  of  the  magnetic  field  for  a  steady-state 
solution.  This  is  exactly  what  happens  experimentally:  strong  oscillations  appear  when 
the  magnetic  field  becomes  too  high;  see,  for  instance.  Fig.  3  of  Ref.[5].  We  come  back 
to  the  limit  rjis  —  0  below. 

8.  A  parameter  difficult  to  choose  is  the  position  of  the  virtual  cathode,  Lep-  Figure  3.8 
shows  that  the  position  of  point  P  affects  more  the  thrust  and  the  efficiency,  when  Epp 
becomes  small.  The  two  efficiencies  involved  in  the  plume,  rjpt  and  rjpc,  behave  in  an 
opposite  way  with  Lpp,  rjpt  decreasing  with  Lpp.  Although  this  issue  requires  further 
analysis,  it  would  seem  desirable  to  choose  Lpp  in  the  range  where  it  affects  the  least  the 
thruster  performances. 

9.  Figure  3.7  shows  that  ap  affects  strongly  the  reverse  ion  flow.  If  we  want  to  keep  a  low 
value  for  that  flow,  there  is  no  much  tolerance  in  the  value  of  unless  other  parameters 
are  changed  simultaneously. 

The  evolution  of  rj  and  rjis  with  and  L,  seems  to  indicate  that  optimum  operation 
conditions  correspond  to  very  low  reverse  flows.  But  operating  the  thruster  in  these  conditions 
has  the  drawback  of  that  a  small  change  in  any  control  parameter  can  lead  the  thruster  response 
into  the  non-stationary  regime.  This  scenario  seems  to  be  the  one  taking  place  in  practice:  the 
thruster  is  operated  near  the  limit  of  the  stationary  regime  and  there  appear  periods  of  a 
strong  oscillatory  response.  If  this  is  the  case,  the  observed  oscillatory  response  would  not 
be  due  to  a  plasma  instability.  For  given  B{x)  and  L,  stationary  solutions  fail  to  exist  for 
low  discharge  voltages  because  the  plasma  temperature  is  too  low  to  produce  enough  ions  to 
maintain  a  supersonic  region  downstream  of  the  ionization  layer  and  quasineutrality  upstream 
of  it.  Longer  channels  and  higher  magnetic  fields  mean  larger  plasma  pressures  and  steeper 
profiles  which  can  be  sustained  steadily  only  by  higher  voltage  discharges. 
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3.7  Discussion 

The  preceding  results  show  that  the  basic  axial  1-D  model  proposed  in  this  chapter  reproduces 
correctly  the  whole  structure  of  the  plasma  in  a  Hall  thruster:  diffusion  region,  ionization 
layer,  acceleration  region,  and  near-plume.  Although  the  main  objective  of  the  model  was  to 
understand  the  plasma  behavior  in  the  thruster,  quantitative  results  for  output  parameters 
like  thrust,  efficiency,  and  discharge  current  agree  well  with  experimental  cases,  which  adds 
reliability  to  the  model. 

Attention  has  been  paid  to  identifying  the  main  phenomena  that  control  the  plasma  re¬ 
sponse.  We  summarize  the  main  aspects  here,  since  previous  models  omit  one  or  another. 
First,  there  is  the  dependence  of  the  ionization  rate  on  the  electron  temperature:  the  strong 
gradient  of  Ri{Te)  around  the  ionization  energy  separates  rather  sharply  the  ionization  and  ac¬ 
celeration  regions  from  the  diffusion  region  and  leaves  this  last  region  with  a  low  temperature. 
Second,  we  have  shown  that  a  reverse  ion  flow  is  forced  by  physical  conditions  around  the  anode: 
solutions  with  v^i  >  0  at  the  anode  boundary  are  not  compatible  with  an  ion-attracting  sheath 
and  therefore  do  not  satisfy  the  condition  of  current  continuity  there.  Since  the  ionization  layer 
starts  around  v^i  =  0,  the  slow  reverse  ion  flow  is  necessary  to  sustain  the  long  diffusion  region, 
also.  The  third  and  main  aspect  to  emphasize  is  the  crucial  role  of  the  electron  pressure  in 
all  plasma  regions:  i)  the  pressure  term  explains  the  subsonic  to  supersonic  transition  of  the 
ion  beam;  ii)  the  electron  pressure  is  the  only  driving  force  for  the  electrons  in  the  diffusion 
region,  where  the  electric  field  is  practically  zero;  iii)  the  pressure  drop,  instead  of  the  diffusion 
losses,  balances  the  gain  in  ion  momentum  in  the  ionization  layer;  iv)  the  maximum  pressure 
is  a  central  parameter  controlling  the  position  of  the  ionization  layer. 

Boeuf  and  Garrigues[ll],  in  spite  of  neglecting  the  electron  pressure,  recovers  a  zero-electric 
field  upstream  the  ionization  region,  by  making  the  B-field  (and,  therefore,  the  diffusion  fre¬ 
quency)  negligible  there.  From  our  preceding  conclusions  we  find  their  solution  debatable  on  the 
following  grounds:  i)  such  a  small  B-field,  which  does  not  correspond  to  usual  settings,  violates 
the  closed-drift  condition  sustaining  the  diffusive  model  of  the  electrons;  ii)  they  neglect  the 
pressure  term  in  the  momentum  equation  in  spite  of  being  much  larger  than  the  electrostatic 
force;  and  iii)  the  compression  work  is  kept,  however,  in  the  energy  equation. 

The  existence  of  a  regular  sonic  point  inside  the  thruster  is  not  just  a  mathematical  artifice. 
Physically,  it  implies  the  local  balance  between  Ohmic  heat  dissipation  into  the  electron  gas  and 
heat  removal  from  it  by  ionization.  The  latter  effect  is  strongly  dominant  near  the  ionization 
peak,  whereas  dissipation,  although  weaker,  dominates  outside  the  layer.  Hence,  the  sonic  bal¬ 
ance  can  only  be  struck  at  a  point  where  ionization  has  nearly  ended  or  has  just  barely  started. 
The  second  option  is  precluded  by  the  fact  that  there  is  no  mechanism  ahead  of  the  ionization 
layer  that  would  provide  the  potential  drop  to  accelerate  the  ions  to  their  forward  speed  of 
sound.  There  is  an  interesting  analogy  between  this  ionization  layer  and  the  shock/combustion 
layer  in  a  strong  detonation  front.  In  both  cases,  the  downstream  exit  point  of  the  layer  is  found 
to  be  sonic,  although,  of  course,  the  layer  inlet  is  supersonic  in  a  detonation  and  subsonic  in  the 
ionization  layer.  The  underlying  mechanism  is  the  same  in  both  problems,  the  condition  of  a 
smooth  sonic  passage  here  being  equivalent  to  the  Chapman-Jouguet  condition  in  combustion 
fronts,  which  states  the  local  balance  between  the  chemical  heating  and  the  heat  diffusion.  One 
interesting  conclusion,  valid  for  both,  detonation  and  ionization  layers,  is  that  the  occurrence  of 
this  ’Chapman-Jouguet  condition’  is  a  consequence  of  the  smallness  of  the  ratio  of  two  kinetic 
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rates,  and  will  only  be  as  precise  as  this  ratio  is  small. 

A  relevant  conclusion,  not  fully  investigated  yet,  is  that  the  presence  of  a  reverse  ion  flow 
seems  essential  for  the  existence  of  stationary  solutions,  the  zero-reverse  flow  condition  indi¬ 
cating  the  transition  to  some  non-stationary  type  of  solutions.  The  parametric  investigation 
on  the  effects  of  the  main  design  parameters  (discharge  voltage,  magnetic  field,  channel  length, 
mass  flow,...)  on  the  thruster  performance  indicates,  for  instance,  that  the  excessive  magnetic 
field  eventually  suppress  the  reverse  ion  flow  and  make  it  impossible  to  maintain  steady-state 
quasineutrality  in  the  diffusion  region,  in  agreement  with  testing  experience. 

Another  issue  requiring  further  research  is  whether  the  transition  from  the  complete  to  the 
incomplete  ionization  modes[17]  corresponds  to  the  transition  from  the  supersonic-exit  to  the 
choked-exit  solutions,  or,  on  the  contrary,  from  the  high-  to  the  low-ionization  branches  of 
supersonic  exits. 

The  main  flaw  of  the  present  model  is  the  extremely  high  peak  temperature.  We  will  try  to 
improve  this  issue  in  Chapters  5  and  6  by  including  plasma  interaction  with  lateral  walls  and 
heat  conduction,  respectively  into  the  electron  energy  balance. 
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Figure  3.1:  [gcerotit]  Regions  leading  to  different  types  of  solutions  between  points  S  and 
B  for  two  values  of  and  a  constant  B-profile.  Line  a  corresponds  to  Gg  =  0;  line  b  to 
ViB  — t  0;  lines  c  and  d  to  Gg  =  0  for  d/  tan  5  =  0.3  and  0.1,  respectively.  Region  2  is  between 
lines  a  and  b;  region  4  (supersonic-exit  regime)  corresponds  to  the  thick  parts  of  line  a,  to 
the  right  of  the  intersections  (asterisks)  with  line  b.  Some  parameter  values  are:  F^,  =  F^, 
=  Sd.O-y/mj/me,  VnA  =  0.1,  CTj  =  2.5,  and  as  =  1/80. 
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Figure  3.2:  [basetit]  Supersonic-exit  solution  for:  xenon,  B{x)  as  shown  in  (a)  \B^  ~  233G, 
Xjn  —  20.3mm,  ~  20.2mm],  the  channel  geometry  shown  in  (b)[L  ~  25mm,  he  —  17.2mm, 

Ac  ~  45cm^,  Lep  25mm],  m  ~  5.32  mg/s,  Vd  ~  293.3  V,  Tgp  ~  2.7  eV,  Vn  ~  300  m/s,  and 
ap  —  1/80.  The  profile  of  the  plume  in  (b)  is  part  of  the  solution;  subplot  (f)  is  a  magnified 
copy  of  subplot(c)  for  region  AS;  asterisks  represent  points  S,  E  and  P;  and  I'd  =  {mc/mi)i>d  in 
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Figure  3.3:  Influence  of  the  discharge  potential  on  the  plasma  response.  Other  parameters  are 
as  in  Fig.  3.2. 
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Figure  3.4:  Influence  of  the  gas  flow  on  the  plasma  response.  Other  parameters  are  as  in  Fig. 
3.2. 
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Figure  3.5:  Influence  of  the  maximum  B-fleld  on  the  plasma  response.  Other  parameters  are 
as  in  Fig.  3.2. 
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Figure  3.6:  Influence  of  the  channel  length  on  the  plasma  response.  Other  parameters  are  as 
in  Fig.  3.2. 
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Figure  3.7:  Influence  of  the  Bohm  diffusion  parameter  on  the  plasma  response.  [Please,  notice 
that  as  in  the  plots  is  as  ■s/’mi/rrif,  in  the  text.]  Other  parameters  are  as  in  Fig.  3.2. 
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Figure  3.8:  Influence  of  the  cathode  distance  on  the  plasma  response.  Other  parameters  are  as 
in  Fig.  3.2. 


Chapter  4 

THE  RADIAL  MODEL 


4.1  Introduction 

The  walls  of  the  discharge  channel  of  the  SPT  class  of  Hall  thrusters  is  made  of  a  ceramic 
insulator.  The  interaction  of  the  plasma  with  these  dielectric  walls  seems  to  affect  greatly 
the  plasma  discharge  and  the  thruster  characteristics[28,  29,  12,  30].  From  simple  reasoning 
several  phenomena  are  expected.  First,  the  zero  current  condition  at  the  dielectric  wall  im¬ 
plies,  in  general,  the  formation  of  an  electron-repelling  sheath  there.  Second,  ions  attracted 
to  the  walls  are  recombined,  thus  reducing  the  actual  propellant  utilization  and  the  thruster 
efficiency.  Third,  energetic  electrons  arriving  at  the  wall  deposit  their  energy  there,  affecting 
the  average  temperature  of  the  plasma  in  the  channel  and,  consequently,  the  ionization  pro¬ 
cess.  Four,  ceramic  materials  present  large  yields  for  secondary  electron  emission.  On  the 
one  hand,  these  electrons  modify  the  sheath  structure,  which  could  attain  charge  saturation 
conditions[31,  32].  On  the  other  hand,  since  these  electrons  are  replacing  primary  electrons, 
which  were  magnetically-guided,  an  extra  near-wall  conductivity  for  the  average  electron  pop¬ 
ulation  has  been  suggested[33]. 

4.2  Presheath/sheath  model 

Here  we  analyze  the  radial  structure  of  the  plasma  in  a  fixed  axial  position  of  the  thruster. 
Fife,  Martfnez-Sanchez,  and  Szabo[29]  were  the  first  to  obtain  consistent  expressions  for  the 
particle,  electron  momentum,  and  electron  energy  wall-source  terms  by  studying  the  behavior 
of  the  space-charge  sheaths  around  the  lateral  walls;  the  effect  of  the  potential  drop  in  the 
radial  presheath  was  added  later  by  Ahedo  and  Martfnez-Sanchez[12]. 

More  recently,  Stephens  and  Ordonez  [31]  and  Jolivet  and  Roussel[32]  have  shown  that  large 
enough  secondary  emission  yields  these  radial  sheaths  can  charge-saturate  the  radial  sheaths, 
a  situation  not  considered  in  the  previous  models  and  with  important  consequences.  Both 
Refs.  [31]  and  [32]  (i)  use  a  particle  approach  (either  kinetic  equations  or  particle-in-cell  for¬ 
mulations),  and  (ii)  force  thermalization  of  the  secondary  emission  at  an  artificial  ’edge  plasma 
source’,  placed  at  the  center  of  the  channel  cross-section.  Our  model  adopts  the  two-scale  (i.e. 
asymptotic),  macroscopic  framework  already  used  in  Ref. [12],  and  includes  trapped  secondary 
electrons  as  a  part  of  the  confined  electron  population  in  the  bulk  of  the  plasma.  In  addition. 
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our  model  includes  diffusion  effects  on  the  ion  radial  motion,  due  to  collisional  processes  and  the 
axial  ion  transport.  The  validity  of  the  macroscopic  approach  for  a  similar  problem  was  shown 
in  Ref.  [34].  With  respect  to  kinetic  models,  a  macroscopic  model  presents  the  advantages 
of:  (i)  following  directly  the  evolution  of  macroscopic  variables  (like  density,  velocity,  ...),  and 
(ii)  being  more  tractable  since  it  solves  ordinary  differential  equations,  exclusively,  instead  of 
integro-differential  ones. 

Figure  4.1  sketches  the  radial  model.  Macroscopic  variables  are  functions  of  r,  and  the  axial 
position  X  appears  just  as  a  parameter  in  the  equations.  To  highlight  the  fundamental  aspects  of 
the  solution,  cylindrical  effects  are  neglected  in  the  present  version;  thus  the  radial  structure  is 
symmetrical  with  respect  to  the  middle-plane  (point  M  in  Fig.  4.1).  In  principle,  the  potential 
profile  is  assumed  monotonic  (in  each  half-channel);  the  transition  to  a  charge-saturated  regime, 
with  local  minima  of  the  potential,  will  be  investigated  afterwards. 

The  model  is  supported  on  three  main  hypotheses: 

(1)  There  are  two  disparate  scales  on  the  radial  profile  of  the  electric  potential:  the  Debye 
length  A(j,  and  the  channel  radial  width  d,  with 

XdCd. 

As  it  is  well  known  from  the  seminal  work  of  Tonks  and  Langmuir [35],  this  justifies  a  two-scale 
analysis,  where,  except  for  the  two  thin  non-neutral  sheaths  tied  to  the  lateral  walls,  the  bulk 
of  the  plasma  is  quasineutral  {the  radial  presheath).  The  electron-repelling  sheaths  are  needed 
to  satisfy  the  zero-current  condition  at  the  walls. 

(2)  Secondary  electron  emission  is  due  mainly  to  the  diffusive  radial  flux,  Pp,  of  primary 
electrons  (population  p)  reaching  the  walls.  As  noticed  in  Ref.  [32],  the  secondary  electrons  can 
be  classified  into  those  becoming  trapped  in  the  channel  and  eventually  thermalized  (population 
s),  and  those  crossing  the  presheath  (in  both  ways)  and  being  collected  back  at  the  opposite 
wall  (two  free  populations  /).  The  present  study  is  restricted  to  the  total  trapping  case,  with 
no  free  populations. 

(3)  A  consistent  treatment  of  the  s-population  is  the  major  challenge  of  the  model.  Sec¬ 
ondary  electrons  flow  continuously  from  the  two  walls  and  are  expected  to  thermalize  through 
multiple  reflections  at  both  sheaths.  Therefore,  in  the  bulk  of  each  cross-section  (i.e.  the 
presheath,  basically)  there  are  trapped  electrons  with  different  time  histories,  and  any  station¬ 
ary  model  must  figure  out  which  is  the  mean  collective  distribution  in  time  of  this  process.  Our 
proposed  model  consists  of:  (a)  treating  primary(p)  and  secondary(s)  electrons  as  independent 
populations  within  each  sheath;  (b)  to  group  all  electrons  into  a  single  electron  population  (e)  in 
the  presheath.  The  details  of  the  trapping  process  are  out  of  the  scope  of  a  macroscopic  model. 
Density,  temperature,  and  diffusion  flux  of  the  population  e  are  defined  through  conservation 
conditions  at  the  two  presheath/sheaths  transition.  We  find  this  model  more  natural  than  the 
one  of  Refs.  [31]  and  [32],  where  secondary  electrons  move  freely  from  the  wall  to  the  center 
of  the  cross-section,  where  they  are  automatically  transformed  into  trapped  electrons  for  the 
other  half-channel.  This  sudden  thermalization  is  rather  artificial  and  introduces  an  undesirable 
electric  field  around  the  cross-section  center  (see  Fig.  1  of  Ref.  [32]).  Another  advantage  of  our 
radial  model  is  to  be  compatible  directly  with  axial  models,  where  a  single  electron  species  is 
generally  considered. 

When  solving  the  equations  for  the  presheath  and  sheath,  we  will  make  use  of  the  funda¬ 
mental  properties  of  this  classical  problem,  like 
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(a)  The  two-scale  analysis  implies  that  the  sheath  is  a  layer  discontinuity  in  the  quasineutral 
scale  of  the  presheath.  The  presheath  has  its  own  spatial  structure,  but  with  spatial  gradients 
infinitely  smaller  than  in  the  sheath.  In  particular,  the  ratio  between  the  electric  fields  in 
presheath  and  sheath  is  A^/d  <C  1. 

(b)  On  the  presheath  scale,  the  transition  to  the  sheath  corresponds,  generally,  to  a  turning 
point  of  the  presheath  solution  (point  Q  in  Fig.  4.2),  that  is,  a  point  where  derivatives  of 
plasma  variables  become  infinity;  in  particular,  the  electric  field  is  infinite  at  point  Q  on  the 
quasineutral  scale.  Indeed,  the  large  gradients  near  the  turning  point  are  announcing  the 
transition  to  the  sheath  with  a  much  thinner  spatial  scale,  the  Debye  length.  The  presheath 
turning  point  is  identified  by  a  certain  condition  on  the  plasma  flow,  which  is  known  as  the 
Bohm  sonic  condition. 

(c)  On  the  sheath  non-neutral  scale,  the  transition  to  the  presheath  corresponds  to  a  zero 
electric  field.  In  general,  the  sheath  is  semi-infinite  in  its  own  scale  and  the  zero  electric  field 
is  attained  asymptotically. 

4.2.1  The  presheath 

Equations  are  written  for  the  half-channel  tm  <  r  <  rw',  the  response  in  the  other  half-channel 
is  symmetric.  We  assume  the  presheath  to  be  constituted  by  two  species  only:  free  ions  (i)  and 
confined  electrons(e).  By  definition,  the  presheath  is  quasineutral  at  every  point: 

rii  ~  rig. 

In  this  issue  we  disagree  with  the  following  two  statements  of  Stephens  and  Ordonez[31]:  (i) 
quasineutrality  is  satisfied  only  globally  in  the  presheath;  and  (ii)  local  quasineutrality  at  the 
presheath/sheath  interface  is  the  matching  condition  there  (when  the  correct  one  is  the  Bohm 
sonic  condition). 

Based  in  2D  planar  equations  for  ions,  we  define  the  following  ID  radial  model: 


d 

(4,1) 

—  TljV>pi 

dr 

dVri  d(j) 

dr  dr 

(4,2) 

where 

1  d 

- (4.3) 

is  the  frequency  for  net  ion  production  in  a  radial  section,  which  includes  both  ionization  and 
the  differential  increment  of  the  axial  ion  flux;  and 

_  I  I  '^xi  ( A  A\ 

Ur  T  I'cx  4  ^  (^•^) 

measures  the  ’friction’  effect  of  collisional  phenomena  and  axial  transport  on  the  radial  velocity 
of  the  ions. 
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Ion  equations  (4.1)- (4.4)  are  completed  with  the  electron  equations 

Vrp  —  Vr 


'-'re  —  ^rt 


dUe  d(j) 

0  ~  -Te-r-  +  erie— . 
dr  dr 


(4.5) 


The  first  equation  results  from  the  zero-current  condition  at  the  lateral  wall.  In  turn,  this  jus¬ 
tifies  neglecting  the  term  VemeneVre  in  the  second  equation  and  Coulomb  collisions  in  Eq.(4.2). 
Also,  the  radial  gradient  of  Tg  (due  to  thermalization  of  secondary  electrons,  for  instance)  has 
been  disregarded. 

Operating  with  Eqs.  (4.1),  (4.2),  and  (4.5),  the  radial  velocity  satisfies 


(Tg  -  rriiV^i) 


dVri 

dr 


(4,6) 


The  transition  to  the  sheath,  point  Q,  corresponds  to  the  turning  point  of  the  presheath  profiles. 
From  Eq.  (4.6),  the  transition  (i.e.  Bohm)  condition  is 


Vri  Vj-jQ  — 


(4,7) 


To  solve  Eq.(4.6),  the  radial  profiles  of  both  frequencies  must  be  specified.  Here  we  will 
consider  the  plausible  case 

—  const,  ^'r(^)  =  const. 

Then,  the  integration  of  Eq.(4.6)  yields 

1  -|-  cr 


arctan(M-v/a)  — 


Vr 


u 


-{r  -  tm), 


(4.8) 


^/^  '^riQ 

where  we  called 

U  Vri  /  '^riQ :  CT  V^  j  Vyj . 

The  profiles  for  the  density  and  the  electric  potential  are  then  determined  from  Eqs.  (4.1), 
(4.5)  and  (4.8), 


e(j)M  —  ecj)  riiM  1  +  i  ,  2\ 

- Tf, - =  In - =  —  ln(l  +  au  ). 

ig  Ui  2a 


(4.9) 


Imposing  the  Bohm  condition  (4.7),  i.e.  u  —  1^  the  potential  drop  (pQM  =  4^m  —  (pq,  and  the 
density  change  across  the  presheath  are  found  to  satisfy 


^4>qm  1  niM  1  +  a  ,  . 

=  In - =  — —  ln(l  +  a). 

ig  riiQ  2a 

Also,  setting  m  =  1  in  Eq.(4.8)  at  rq  —  ru  —  tw  —  fM  —  d/2,  we  have 

Vr  A 

—  =  Vr[a)  —  2 1  — ^  arctan  Ca  —  1  , 

Vq  ^  y/ a  ) 

Vr{(x) 


(4.10) 


(4,11) 


-  =  = 

Vo  a 
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which  is  the  parametric  equation  relating  the  two  relevant  frequencies  of  the  radial  problem. 
Here, 


(4.12) 


is  the  frequency  of  reference  for  radial  processes.  Equation  (4.11)  expresses  that  the  produc¬ 
tion/loss  frequency  is  not  free.  The  dielectric  wall  implies  the  presence  of  a  non-neutral  sheath, 
which  forces  the  Bohm  condition  on  the  ion  flow  to  the  wall,  and  this  makes  Vyu  not  free.  Its 
actual  value  depends  on  the  acceleration  process  for  ions  and  this  is  the  meaning  of  the  rela¬ 
tionship  Vwi^r)-  In  the  model  of  Ref.  [12],  where  we  assumed  Vy  —  0,  Vyj  was  totally  determined 
and  equal  to  4/3. 

Figure  4.3  shows  the  profiles  ^(r)  for  three  values  of  Vy.  On  the  presheath  scale,  the  jumps 
of  the  different  plasma  variables  across  the  sheath  are  vertical  discontinuities  starting  at  the 
turning  point  Q,  and  are  to  be  determined  from  the  sheath  equations. 

Figure  4.4  shows  the  evolution  of  Vw  and  e(t)QM/Te  versus  Vr-,  Eqs.  (4.10)  and  (4.11).  Asymp¬ 
totic  expressions  for  small  and  large  values  are 


Ur  <.1  ■. 

4 

7/  ~  - 

“'“3’ 

^4>qm 

T 

e 

1 

“  2’ 

Z>r  >  1  : 

7/  ~  - 

~  5 

^(pQM 

Ur 

~  In  — 

Ur 

Te 

TT 

Physically,  a  large  ’friction’  implies  a  large  potential  drop  to  accelerate  the  ions  to  the  Bohm 
velocity.  As  a  consequence,  is  low,  the  plasma  current  to  the  wall  is  low,  and  i>yj  is 

low. 


4.2.2  The  sheath 

Both  ion  and  electron  dynamics  present  changes  in  the  sheath.  For  Xa/d  0  the  sheath 
thickness  is  negligible  compared  to  the  channel  width,  the  contribution  of  Vyj  and  Vy  to  ion 
dynamics  within  the  sheath  is  negligible,  and  Eqs. (4.1)  and  (4.2)  yield  the  conservation  of  ion 
particles  and  energy.  Within  the  sheath,  electrons  are  classified  into  the  main  primary  popu¬ 
lation  (p),  which  diffuse  towards  the  wall  and  causes  secondary  emission,  and  the  population 
(s)  of  secondary  electrons,  which,  after  crossing  the  sheath,  become  trapped  in  the  presheath, 
in  a  process  not  analyzed  here.  The  set  of  macroscopic  equations  for  the  three  species  in  the 
sheath  are: 

UiVri  =  const  ^  gi  =  riiq 

IrriiV^i  +  e(j)  —  const  =  |Te  -h  e^g, 

UpVrp  —  const  —  Qp  —  npw^jTp/2'Kmf., 

A  A  (4-13) 

e4)  -  e(pQ 
np  =  UpQ  exp - - - , 

J-p 

risVrs  =  const  =  g^, 

—  ecj)  —  const  =  Tyj. 
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The  value  of  the  primary  flux,  Qp^  was  calculated  assuming  a  Maxwellian  distribution  of  tem¬ 
perature  Tp.  The  mean  energy  of  the  population  of  secondary  electrons  when  they  leave  the 
wall  is  Tyu,  which  can  be  assumed  of  the  order  of  the  wall  temperature.  Therefore,  one  expects 

fy,/Tp  <  1, 

which  allows  one  to  work  in  the  asymptotic  limit  Tyj/Tp  0,  except  in  the  very  vicinity  of 
point  W.  Finally,  the  low  value  of  and  the  acceleration  of  the  secondary  electrons  as  they 
cross  the  sheath  justifies  neglecting  thermal  effects  for  them. 

The  flux  of  secondary  electrons  is  assumed  proportional  to  the  flux  of  primary  electrons: 

9s  =  -5yj9p  {9s  <  0),  (4.14) 

with  5yj  the  effective  secondary  emission  yield.  In  addition,  the  zero-current  condition  at  the 
wall  relates  the  electron  and  ion  fluxes. 


9i^  9p  +  9s  ^  {'^-  ^w)9p-  (4.15) 

To  compute  the  effective  yield  for  secondary  emission  we  start  with  the  case  of  a  monoener- 
getic  electron  beam.  If  E  is  the  beam  energy,  we  average  over  the  different  angles  of  incidence, 
and,  for  E  low  enough,  we  consider  that  dyj  satisfies 

^w{E)  ^  (4-16) 

Here,  Eyj  is  the  energy  for  100%  yield,  and  depends  strongly  on  the  wall  material;  the  typical 
range,  for  the  materials  of  interest  in  a  Hall  thruster,  would  be  Eyj  ~  20  —  100  eV.  The  exponent 
p  is  around  0.5,  for  the  same  materials.  Next,  for  a  Maxwellian  population  of  temperature  Tp, 
the  effective  secondary  emission  yield  can  be  computed  by  averaging  the  law  (4.16)  over  the 
different  electron  energies.  This  yields[29] 

sun)  -  {^y.  0-17) 

with 

P  _  Eyj 

r(2  +  p)i/p’ 

for  p  —  0.5,  one  has  Eyj  ~  0.57Eyj.  For  Eyj  ~  30  eV  (Boron  Nitride  ceramic  referenced  in  [12]) 
one  obtains  Eyj  ~  16.64  eV  for  p  —  0.576  ,  or  F'.u,  ~  17.0  eV  for  p  —  0.5.  For  Eyj  ~  70.2  eV 
(BNAIN  ceramic  referenced  in  [32]),  one  has  Eyj  ~  40.0  eV  for  p  —  0.5. 

Plasma  equations  (4.13)  are  completed  with  Poisson  equation 


Co  cPcj) 
e  dr'^ 


^Up  +  ris-  rii. 


Integrating  this  equation  together  with  Eqs.(4.13)  and  imposing  that  the  electric  field  goes  to 
zero  asymptotically  at  point  Q,  that  is 


djeHTe) 
d{r/Xd)  Q 
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one  finds  the  conservation  equation 

j[%)‘  =  V(e^)-V{e<j,Q),  (4,18) 

where 


U{e<j))  =  UpTp  +  m^nsvl^  + 


is  the  Sagdeev’s  potential.  Equation  (4.18)  expresses  the  balance  of  the  dynamic  and  electric 
pressures  of  the  whole  plasma  across  the  sheath.  Notice  that  the  first  derivative  of  f/  measures 
the  local  electric  charge:  —  np  +  Ug  —  rii.  Also,  the  electric  field  at  the  wall  measures  the 

total  electric  charge  in  the  sheath: 


Co  d(p 


w 


(4.19) 


Two  boundary  conditions  at  point  Q  are  related  to  U{e(f)).  First,  plasma  quasineutrality  is 
achieved  there. 


U'q  =  -  n-iQ  —  0. 


(4.20) 


Second,  with  the  presence  of  two  electron  populations  in  the  sheath,  the  Bohm  condition  takes 
the  form 


U: 


n  _  "^pQ 


n 


sQ 


n. 


iQ 


Q 


^eV;gQ 


=  0. 


rriiV, 


(4.21) 


riQ 


Finally,  to  match  correctly  the  solutions  of  presheath  and  sheath,  we  must  define  the  rela¬ 
tionships  among  the  magnitudes  of  the  different  electron  populations  at  both  sides  of  point  Q. 
These  conditions  are 


"neQ  —  npQ  -h  UsQ, 

9eQ  9pQ  T  9sQ  (1  ~  ^w)9pQ^ 
f^eQ  _ 


(4.22) 


nrieV. 


rsQ 


The  last  expression  makes  equivalent  the  two  forms  of  the  Bohm  condition,  Eqs.(4.7)  and  (4.21), 
and  determines  Te(Tp).  One  can  see  that  Tg  >  Tp,  as  expected  from  the  grouping  of  the  primary 
and  secondary  populations  into  a  single  one. 

Notice  that  e—  and  p—  magnitudes  coincide  for  5^  —  0,  as  it  should  be.  Let  us  use  this 
case  to  justify  that  the  conditions  defining  the  transition  presheath/sheath  (i.e.  point  Q)  are 
unique.  On  the  one  hand,  a  solution  for  the  presheath  exists  only  for  v^iQ  <  \fT^Jmi\  on  the 
other  hand,  a  solution  in  the  sheath,  around  point  Q,  requires  that  Uq  >  0,  which  (for  —  0) 
states  that  Vriq  >  ^/Te/rrii.  Therefore,  Vriq  —  ^jTj^i  is  the  only  valid  transition  condition. 

The  sheath  equations  and  boundary  conditions  have  shown  that  Tg,  and  rigg  fully  char¬ 
acterize  the  sheath  solution  (except  for  the  marginal  parameter  f^/Tg  <C  1).  Furthermore, 
using  Tg  and  rigg,  and  the  related  Debye  length 


Aw  — 


CUeq, 
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to  non-dimensionalize  the  sheath  equations,  the  dimensionless  solution  depends  on  just  the 
secondary  emission  yield,  5^.  In  particular,  the  inner  spatial  variable  to  be  used  is 


y  = 


r-rw 


(4.23) 


Taking  into  account  that  Uq,Uq  —  0,  the  potential  around  point  Q  comes  from  the  expansion 


of  Eq.(4.18),  with  Uq  <  0  for  an  ion-attracting  sheath.  For  each  S^,  the  potential  profile  (f){y) 
in  the  sheath  is  obtained  by  determining,  first,  the  sheath  potential  jump,  (pwQ  —  <pQ~<Pw->  and 
performing,  then,  a  numerical  quadrature  of  Eq.(4.18),  from  the  vicinity  of  point  Q  towards 
point  W,  using  the  above  expansion  as  initial  condition. 

Figure  4.5  shows  the  potential  profile  in  the  sheath  (between  points  W  and  Q,  with 
yq  —  — oo)  for  three  values  of  Thus,  the  whole  potential  profile,  from  the  channel  center 
to  the  wall,  depends  on  the  values  of  5^  and  and  is  obtained  by  matching  the  respective 
presheath  and  sheath  profiles. 

Figure  4.6  shows  the  evolution  of  the  main  sheath  parameters  with  5yj  up  to 


~  0.983,  (4.25) 

when  the  electric  field  at  the  wall  boundary  is  zero,  d(j)/dy\Q  —  0,  which  corresponds  to  the 
charge  saturation  limit  (CSL),  to  be  discussed  below.  Fig.  4.6  shows  that  Upg/ueg  and  Tpg/Tgg 
remain  very  near  unity  up  to  ~  0.9,  and  then  the  contribution  of  the  secondary  emission 
grows  to  a  10-15%  for  the  CSL.  This  means  that  the  properties  of  the  electron  population  in 
the  presheath  are  mainly  determined  by  primary  electrons,  which  adds  reliability  to  our  model 
of  a  unique  electron  population  in  the  presheath. 

For  Sya  —  0,  the  sheath  potential  is  large,  e(f)wg/Te  ~  5.28.  For  Sy^  >  0,  the  sheath  potential 
satisfies 

Since  the  last  term  on  the  right  side  is  small  for  any  Syj,  the  first  term  on  the  right  side 
dominates  until  the  secondary  yield  satisfies  1  —  Syj  ^  10^^,  that  is  up  to  the  vicinity  of  the 
charge  saturation  limit.  This  explains  why  the  potential  profile  does  not  change  significantly 
until  Syj  does  not  approach  the  CSL.  The  potential  drop  in  that  limit  is  ~  0.88. 


(4.26) 


4.2.3  The  charge  saturation  regime 

As  the  secondary  emission  yield  increases,  the  electron  density  near  the  wall  increases  and  the 
electric  field  at  the  wall  decreases,  as  shown  in  Fig.  4.6,  until  it  becomes  zero  at  the  CSL.  From 
Eq.(4.18),  the  condition  for  charge  saturation  is 


0  —  U{e(j)w)  —  U{e(j)g), 
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or,  in  terms  of  the  dynamic  pressures  of  the  different  species, 


Tpi^TlpQ  IT'pw')  A  ^T^edsiA'sQ  ^sw')  ^zq')-  (4.27) 

Since  the  electric  field  at  the  wall  depends  just  on  6^,  the  charge  saturation  limit  corresponds 
to  a  certain  value  of  that  parameter,  as  we  saw  before. 

To  explain  the  value  obtained  for  at  the  CSL,  Eq.(4.25),  let  us  recall  the  classical  CSL 
problem  for  planar  sheaths  with  e(j)wQ/Tp  ^  1  and  Qs  independent  of  Qp.  For  them,  Eq.  (4.27) 
simplifies  into  rriigiUiw  —  i^edsUsQ-,  where  from  the  ion  and  electron  counter-streaming  fluxes 
satisfy 


9s 


(4.28) 


which  is  known  as  the  Langmuir  condition.  Physically,  this  law  is  explained  by  two  facts:  on  the 
one  hand,  integrating  Poisson  equation  the  total  electric  charge  in  a  charge-saturated  sheath  is 
zero,  so  the  average  densities  of  accelerated  ions  and  electrons  are  the  same;  on  the  other  hand, 
the  average  kinetic  energies  of  ions  and  electrons  are  the  same  and  of  the  order  of  the  potential 
jump.  For  e(j)wQ/Tp  —  0(1),  more  terms  in  Eq.(4.27)  will  contribute  to  the  fluxes  ratio  but  we 
can  expect  the  Langmuir  condition  (4.28)  to  be  still  verified,  except  for  a  factor  of  order  one. 
In  our  case,  applying  both  the  Langmuir  and  the  zero-current  conditions,  the  CSL  is  expected 
to  happen  for 


ys 

in  agreement  with  the  value  given  by  Eq.(4.25). 

Let  us  discuss  now  the  plasma  response  once  the  CSL  has  been  attained.  For  the  real  case 
Tu,  >  0,  a  secondary  emission  yield  larger  than  leads  to  a  potential  profile  with  a  minimum 
near  the  wall  (at  a  point  W’,  let  us  say),  which  acts  as  a  potential  barrier  to  turn  back  to 
the  wall  the  excess  of  secondary  flux.  The  actual  secondary  flux,  reaching  the  presheath  is 
the  flux  at  point  W’.  Since  d(l)/dy\w'  —  0,  the  secondary  flux  between  W’  and  Q  satisfies  CSL 
conditions,  that  is  Eqs.  (4.14)  and  (4.25) 


9's  =  -^*yj9p- 

The  typical  size  of  the  potential  barrier,  to  turn  back  the  excess  of  flux,  —  ^^1,  is 

e(^iiv  —  (pw)  ~  Tyj. 

Then,  in  the  asymptotic  limit  T^/Tp  — >■  0,  the  region  W’W  is  negligible  and  the  CSL  (in  Fig. 
4.6,  for  instance)  instead  of  identifying  a  single  operation  point,  identifies  a  whole  domain 
of  operation,  the  charge-saturated  regime  (CSR),  satisfying  5yj  —  at  all  operation  points. 
Therefore,  the  unified  expression  of  the  secondary  yield  for  the  non-CSR  and  CSR  is 

5^{Tp)  =  minj^;, 


(4.30) 
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Although  Fig.  4.6  gives  a  complete  information  of  the  sheath  solution  for  the  two  regimes, 
a  clearer  view  of  the  change  of  behavior  at  the  SCL  is  given  by  Fig.  4.7,  which  shows  the 
evolution  of  (pwQ  with  Tg  for  different  E^,;  one  has 

~  0.88 


at  all  points  of  the  CSR. 

A  final  but  very  relevant  consequence  of  the  existence  of  the  charge-saturation  regime  is 
that  the  electron-repelling  sheath  exists  always,  contrary  to  what  we  expected  previously  in 
Ref.  [12].  The  charge-saturated  sheath  will  limit  particle  and  heat  losses  in  the  high  Tg  range. 


4.3  Wall-source  terms  in  the  axial  model 

Once  the  radial  profiles  are  known,  we  can  evaluate  the  source  terms  of  Eqs.(2.21)-(2.24).  The 
quasiplanar  approximation  will  be  used. 


4.3.1  Particle  losses 

The  ion  and  electron  (negative)  gains  per  unit  of  volume  in  a  radial  section  are,  Eq.  (2.21), 


or,  using  Eq.  (4.1), 

2 

Sw  -  -  ^  /  niiy^dr  =  (4-31) 

«  JM 


where  hi  —  hg  is  the  average  plasma  density  in  a  radial  section.  From  the  two  above  expressions 
one  has 


hj  _  2  riiQ 
f^iM  hyj  TliM 

and  hi/riiM  depends  only  on  hr  and  is  of  order  unity. 
The  production  frequency  satisfies 


(4.32) 


(4.33) 


with  [Eq.(4.12)] 


and  Vwi^r)  given  by  Eq. 
in  Eq.  (4.4),  by 


/15mm\  /Tg\i/2  , 

)U)  x2xl0=s-‘ 

(4.11).  The  scaling  law  Vri  oc  -s/Tl,  suggests  to  approximate  hr,  defined 


ly, 


(x)  ~  —(vi 

Vq  \ 


+  ^cx  + 


'^xi  dy/Tg^ 


(4.34) 
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an  expression  that  can  be  made  independent  of  r.  The  last  term  includes  the  ion  Mach  number 
of  the  axial  flow  which  is  large  near  the  exit.  Therefore  the  slope  of  Te{x)  affects  strongly  the 
particle  losses  there. 

The  losses  in  plasma  current  (per  unit  of  axial  length)  can  be  calculated  from 


cS'iijAfi 


f  Ac  \  /  he  \  ^ 

vdOcm^/  vlO^®m^^/ 


lO^s-i 


)  X  0.64A/ 


cm, 


with  Ac  the  cross-section  area  of  the  channel  (the  reference  values  used  here  are  typical  of  a 
SPT-100  thruster).  Unless  the  friction  parameter  is  large,  particle  losses  are  going  to  be 
high. 


4.3.2  Near-wall  conductivity 

The  gain  per  unit  of  volume  in  azimuthal  momentum  of  the  electrons  in  a  radial  section,  due 
to  the  fluxes  of  primary  and  secondary  electrons,  Eq.(2.23),  is 

2  2  ^ 

for  fyj/Tp  -)•  0. 

The  equations  for  electron  current  follow  the  closed  drift  approximation.  Taking  into  account 
radial  effects,  they  can  be  written  as 

0  =  erie- - — UeTe  -h  UJcmcUcVec, 

ax  ax 

0  =  -mcflciuJeVxe  +  ^eVOe)  + 

where  the  azimuthal  ’friction’  force  due  to  the  walls,  satisfies 

^wde  ~  ^w6e  aicSyjc^Se  ^wm'^e'^e^dej  (4.35) 


with 


^wm 


(4.36) 


the  related  frequency,  and 


Syj 

l-Sy, 


an  enhancement  function,  shown  in  Fig.  4.8,  due  to  secondary  emission,  which  can  multiply 
the  production  frequency  by  a  factor  up  to  60  in  the  CSR. 

The  ’friction’  frequency  adds  a  second  term  to  the  electron  axial  current,  which  now  satisfies 


-  -  +  ^wm  (  ^  d,Pc\  ,  , 

tncVx^c  —  — - e [eUcE^  ~  (4.37) 

This  contribution  is  the  ’near- wall  conductivity’ (NWC)  of  Ref.  [33].  Here,  the  NWC  does  not 
present  radial  oscillations  along  the  cross-section,  because  Vyjm  represents  an  average  value. 

Since  typical  values  of  Vc  are  of  order  lO^s^^,  the  near- wall  conductivity  will  be  a  marginal 
effect  except  for  the  charge-saturated  regime  and  a  large  friction  parameter,  i>y. 
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4.3.3  Energy  losses 

The  (negative)  gain  in  electron  internal  energy,  per  unit  of  volume,  due  to  wall  interaction, 
satisfies 

Qwe-—^\[2^eTeVre  +  qre\j^-  Vre—^dvj,  (4.38) 

with  Qj-e  accounting  formally  for  heat  conduction.  The  computation  of  the  different  terms  yields 

Q'we  =  -J^wefleTe,  (4.39) 

where 


^we  —  T  /^e(^y) 


(4.40) 


is  the  related  frequency,  and 


3a  +  l 

v?du 


are  two  enhancement  functions.  The  first  one,  f^eidw),  plotted  in  Fig.  4.8,  depends  on  the 
sheath  characteristics;  it  goes  from  /fg  —  7.28,  for  6^  —  0.,  /dg  —  jd*  ^  104.5,  for  6^  —  The 
second  one,  fd'^{i>r)-,  depends  on  the  presheath  characteristics  and  has  a  lesser  effect  in  we 
find  /3'  <C  /dg  for  Vr  —  0(1),  and  Vwfd'e  ^  i^o(lnf'r)/f'r  ^  for  Vr  1.  For  the  sake  of  clarity 
we  suggest  to  use  the  expression 

1  65 

Pe{dw)  ^  5-62  +  '  ,  (4.41) 

l-dw 


which  approximates  the  exact  function  with  an  error  less  than  10%. 

The  losses  in  internal  energy  (per  unit  of  axial  length)  can  be  calculated  from 


-  (ii^)  (») 


Tg\3/2 


Fig.  4.9  computes  \Q'^^Ac\  in  terms  of  Tg  for  different  values  of  the  characteristic  yield  energy, 
Eyj.  As  in  the  case  of  the  sheath  potential  (pwQ^  Fig-  4.7,  there  is  a  change  in  behavior  from 
the  non-CSR  to  the  CSR.  For  a  given  Tg  and  below  the  charge-saturation  regime,  energy  losses 
increase  with  Ey^  decreasing.  Within  the  charge-saturation  regime  energy  losses  are  independent 
of  Eyj  and  grow  with  Tg^^-  Figure  4.9  shows  that  the  energy  losses  can  be  enormous  if  there 
is  not  a  large  ion  friction,  or  any  other  phenomenon  able  to  reduce  significantly  the  electron 
current  to  the  wall. 

For  wall  heating  considerations,  the  important  magnitude  is  the  total  energy  deposited  by 
electrons,  which  is  the  sum  of  the  heat  loss,  Q'yjg-,  and  the  loss  of  azimuthal  energy,  rrigvlJigVyj. 
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4.4  Conclusions 

We  have  presented  a  model  of  the  radial  structure  of  the  plasma  interaction  with  dielectric  walls, 
supported  in  three  pillars:  (i)  a  two-scale  asymptotic  analysis;  (ii)  a  new  model  to  account  for 
trapped  secondary  electrons;  and  (iii)  the  inclusion  of  ion  radial  ’friction’,  due  to  collisional 
processes  and  axial  transport.  The  presheath  properties  depend  on  the  ion  friction  parameter 
i>r,  whereas  the  sheath  properties  depend  on  the  effective  secondary  emission  yield,  Sw  A  large 
friction  increases  the  presheath  potential  drop,  and  reduces  the  plasma  current  to  the  wall. 
The  charge  saturation  limit  is  attained  for  ~  0.983,  and  marks  the  transition  to  a  charge- 
saturated  regime,  with  a  constant  ratio  for  (sheath  potential  drop)/ (plasma  temperature).  As 
a  consequence,  the  electron-repelling  sheath  never  disappears. 

Particle  exchanges  at  lateral  walls  lead  to  losses  in  particles,  energy,  and  azimuthal  momen¬ 
tum.  The  near- wall  conductivity  has  been  determined  from  the  loss  in  azimuthal  momentum.  It 
seems  to  be  a  moderate  effect  except  for  operation  near  or  under  the  charge-saturation  regime. 
For  typical  thruster  conditions  the  model  predicts  too  large  losses  of  particles  and  energy,  un¬ 
less  a  large  ion  ’radial  friction’  is  assumed.  The  coupling  of  the  radial  and  axial  models,  now 
underway,  will  show  whether  a  large  radial  friction  is  attainable. 

Total  trapping  of  secondary  electrons  does  not  seem  the  case  in  practice,  according  to  typical 
values  of  the  total  electron  collision  frequency.  Partial  trapping  implies  a  double  population 
of  secondary  electrons  traveling  freely  from  wall  to  wall.  This  more  general  case,  now  under 
study,  will  certainly  modify  the  radial  structure  and  could  provide  and  additional  mechanism 
to  inhibit  radial  transport  to  the  walls. 


CHAPTER  4.  THE  RADIAL  MODEL 


52 


Figure  4.1:  Sketch  of  the  radial  model  in  the  quasiplanar  case.  Point  M  represents  the  mid-plane 
of  the  channel.  The  sheath  is  the  region  of  large  gradients  of  ^(r). 


Q  W 

Figure  4.2:  Two-scale  analysis:  point  Q  is  the  presheath/sheath  transition  with  d<f)/dr  — oo  in 
the  quasineutral  scale,  gp  and  Qs  are  the  fluxes  of  primary  and  secondary  electrons  at  the  wall; 
in  the  presheath,  both  populations  are  grouped  in  a  single  one  with  ge  —  9p  +  gs- 
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Figure  4.3:  Presheath:  Potential  profile  for  different  values  of  the  dimensionless  friction  param¬ 
eter. 
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Figure  4.4:  Presheath:  Influence  of  the  radial  friction  frequency  on  the  production  frequency 
and  the  potential  drop  across  the  presheath. 


Figure  4.5:  Sheath:  Potential  profile  for  different  values  of  the  secondary  emission  yield.  For 
—  0-983  it  is  d(j)/dy\w  —  0,  marking  the  charge-saturation  limit.  Point  Q  is  here  at 
Vq  =  -oo- 
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Figure  4.6:  Sheath:  Influence  of  the  secondary  emission  yield  on  potential  drop,  electric  held  at 
the  wall,  and  electron  temperatures  and  densities.  Point  CSR  represents  the  charge-saturation 
regime;  e(j)wQ/Te\csR  -  0.879,  Tp/Te\csR  -  0.865. 
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Figure  4.7:  Evolution  of  the  sheath  potential  drop  with  the  plasma  temperature  for  different 
ceramic  materials  (i.e.  E^).  The  straight  line  corresponds  to  the  CSR. 


Figure  4.8:  Enhancement  functions  for  the  effective  frequencies  accounting  for  near-wall- 
conductivity,  PwrYi,  and  energy  losses,  lywe- 


Figure  4.9:  Energy  losses  per  unit  of  channel  length  for  different  wall  materials,  Ac  —  40  cm^, 
d  —  15mm,  hg  =  and  —  1-  The  solid  line  corresponds  to  the  CSR.  Energy  losses 

are  proportional  to  Pw 


Chapter  5 

EFFECTS  OF  LATERAL  LOSSES 


5.1  Governing  equations 


Here  we  add  the  effects  of  the  plasma  interaction  with  the  thruster  walls,  studied  in  the  previous 
chapter,  into  the  basic  steady-state  model  treated  in  Ch.  3.  For  sake  of  clarity  in  the  exposition 
and  to  avoid  continuous  references  to  previous  chapters,  we  give  a  brief  but  complete  account 
of  the  model  equations. 

The  stationary,  macroscopic  equations  for  the  quasineutral  plasma  between  the  entrance  to 
the  anode  sheath  (point  B)  and  the  external  neutralization  surface  (point  P)  are 


\-^{AneV^i)  =  \-^{AneV^e)  =  -\-^{AnnVn)  =  neiPi  -  Uyj), 

A  ax  A  ax  A  ax 

^  ^  /  A  2  \  \ 

^w'^xih 

d  d(p 

0  =  -  —UeTe  +  eUe- - UdmeUeV^e, 

ax  ax 

~A~dx^  ^d^e^e'^xe  U-ujgTlgTe 


The  ionization  frequency  follows 

^rjRi(Te), 

the  axial  diffusion  frequency  for  the  magnetized  electrons  satisfies 


OJZ 


V.  -h  Pu 


with 


(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 


(5.6) 


—  ^en  +  ^ei  + 

the  electron  collision  frequency,  grouping  contributions  from  e  —  n  and  e  —  i  collisions,  and  Bohm 
anomalous  diffusion;  ^wm  and  Py,f.  are  frequencies  accounting  for  particle  losses,  near- wall 
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conductivity,  and  energy  losses  at  the  lateral  walls;  and  Vnw  is  the  axial  velocity  of  ions  after 
wall  recombination  and  accommodation. 

The  area  A{x)  and  thickness  d{x)  of  the  plasma  jet  cylindrical  cross-section,  are  constant 
within  the  channel  and  satisfy 


—  =  -  tan  5,  X  >  L,  (5.7) 

ax  a 

in  the  plume,  with  5  the  local  mean  angle  of  divergence  of  each  boundary,  to  be  defined  below 
and  a  =  In  ^  const  =  In  d  -|-  const. 

The  azimuthal  components  of  the  electron  velocity  and  the  heat  flux  follow 

'^de  Qde 

_  _ 

'^xe  Qxe 


(5.8) 


with  (ne/z^e  ^  1- 


5.1.1  Frequencies  for  wall-source  terms 

We  summarize  here  the  expressions  obtained  in  Ch.4.  To  help  solving  the  axial  model,  some 
approximations  are  made  too.  The  particle  loss  frequency  satisfies 


— 3 - 


(5.9) 


with  vq  the  reference  frequency  for  the  radial  motion,  and  Vyj  a  dimensionless  factor,  which 
depends  on  the  characteristics  of  the  radial  motion.  In  this  work  we  will  take  the  simplest  case 


Vyj  —  const. 


Future  work  will  have  to  take  into  account  the  influence  of  the  radial  dynamics  of  ions  and 
secondary  electrons  on  that  is  the  functional  relation  discussed  in  Ch.4. 

The  frequencies  for  near-wall  diffusion  and  heat  losses  are 


^wm  Pmi.^w') 


1  -  Sy,  ’ 


^we  —  ^  5.62  -|- 


1.65 

1-Sy 


(5.10) 


with  j3yn  and  j3e  enhancement  functions,  which  depend  on  the  effective  secondary  emission  yield 
dw{Te).  Assuming  that  the  secondary  emission  yield  for  an  impinging  monoenergetic  electron 
beam  of  energy  E,  follows  the  law 


5yj{E)  ~  ^JE/Eyj, 


the  effective  yield  for  the  quasi-Maxwellian  population  of  primary  electrons,  satisfies  Eq.(4.30) 


5^(Te)  ~min|5;,y^}. 
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with 


El, 


the  temperature  leading  to  100%  secondary  emission[36]  and 


5;  ~  0.983 

the  upper-bound  value  of  5^,  corresponding  to  the  charge-saturated  regime  for  the  lateral 
sheaths.  Therefore,  the  maximum  values  of  the  enhancement  functions  are 


/3;  ^  60,  p:  ^  105, 


and  the  charge-saturation  regime  corresponds  to  Tg  >  0.967£'.u,. 


5.1.2  Singular  points  and  boundary  conditions 

Solving  Eqs.  (5.1)-(5.5)  for  the  spatial  derivatives  we  obtain  a  matrix  relation  of  the  form 

dY 

(1-M^)— =  F(r),  (5.11) 

where  Y  —  (Tg,  rig,  v^i,  v^e,  Vn,  (f)  groups  the  7  plasma  variables,  F  is  a  regular  function  and 

M  =  (5.12) 

y/5Te!3mi 

is  the  Mach  number  for  the  ion  axial  flow.  For  instance,  the  equations  for  rig  and  v^i  can  be 
written  as 


1  drig 
rig  dx 

with  P  —  Te{l  —  M^)  and 


G  dVxi  / G  da\ 


(5.13) 


^xeO 


r2  3 

—  Oi-iE-i  Tg  —  771^1! xe{‘^'Vxi  '^n) 


2  3  2  d  In  ^ 


G, 


(5.14) 


Sonic  points  and  boundary  conditions  are  discussed  in  Chap. 3.  The  seven  boundary  condi¬ 
tions  needed  here  are: 

i)-ii)  The  injected  flow  of  neutrals  at  the  anode,  m,  and  their  velocity,  VnB  —  VnA,  are  known. 

iii)  The  electron  temperature  at  the  neutralization  surface,  Tgp^  is  known. 

iv)  The  potential  drop  between  points  A  and  P  is  the  discharge  voltage,  (I)a  —  (f^p  —  Va- 

v)  The  presheath/sheath  transition  at  point  B  requires  Mb  —  —1,  which  is  the  Bohm  sonic 
condition  on  the  ion  back-flow. 

vi)  The  potential  jump  at  the  anode  sheath  satisfies  Eq.(3.11) 


^4>AB  1  CeB 

- 

AeB  4  VxeB 
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with  Ce  =  y/STefTVme- 

vii)  There  is  a  regular  sonic  point  inside  the  channel  (point  S),  which,  according  to  Eq.(5.13), 
is  characterized  by 


Gs  —  0  at  ALg  —  1. 


(5.15) 


In  addition,  we  must  define  an  expression  for  the  angle  of  divergence  6.  Following  Ref.  [37] 
we  take 


y^bTeElArrii 
tan  0  —  — - 

'^xi 

Notice  that  discharge  current.  Id,  position  of  forward  sonic  point,  xs,  and  divergence  angle 
at  the  channel  exit,  Se,  are  part  of  the  solution. 

5.2  Axial  plasma  structure 

Figure  5.1  shows  the  first  complete  results  we  have  obtained  with  the  preceding  model  of  radial 
losses.  Profiles  of  most  plasma  variables  are  shown  for  a  SPT-100  type  of  thruster.  The  solution 
shown  here  presents  moderate  wall  losses.  However,  the  values  selected  for  the  parameters  that 
determine  radial  losses:  ~  0.17  and  Ey^  —  lOOeP,  were  aimed  mainly  to  test  the  possibility  of 

obtaining  valid  solutions.  Our  procedure  to  obtain  these  solutions  is  to  carry  out  a  parametric 
continuation  from  the  solution  without  wall  losses  Uy,  —  0.  Actual  values  of  Eyj  for  BN-based 
materials  are  in  the  range  Eyj  ~  2Ey,jZ  ~  15-40  eV.  From  Eq.(4.11)  Vy,  ~  0.17  corresponds  to 
Vy  ~  50,  and  ~  lO^Hz,  which,  from  the  definition  of  Eq.(2.59),  seems  too  large;  however, 
partial  trapping,  not  investigated  yet  could  reduce  Vy,.  Attempts  to  increase  Vy,  and  reduce  Ey,, 
beyond  the  case  of  Fig.  5.1,  have  encountered  different  convergence  difficulties.  One  of  them  is 
the  appearance  of  complex  values  for  the  plasma  derivatives  at  point  S  preventing  the  launch 
of  the  integration. 

The  solution  of  Fig.  5.1  shows  that  energy  losses  at  lateral  walls  smooth  the  temperature 
profile,  the  maximum  temperature  decreasing  from  about  90  eV  with  no  wall  losses,  to  about 
68  eV  now,  which  corresponds  to  a  maximum  secondary  yield  of  by,  ~  0.82,  far  (in  terms  of 
losses)  from  the  charge-saturation  value  b*^^.  Energy  losses  take  place  in  the  acceleration  region, 
where  Tg  is  high.  On  the  contrary,  losses  in  plasma  current  are  concentrated  in  the  rear  part 
of  the  channel.  The  (relative)  ion  back-flow  at  the  inner  boundary  of  the  ionization  layer  is 
rjiE  —  —0.37,  whereas  at  the  anode  is  just  rjis  —  —0.05,  which  means  a  current  loss  in  the 
diffusion  region  of  1.2  A.  The  total  plasma  current  deposited  at  lateral  walls,  f  Sy,dx,  with 

Syj  CTl^l^yjAQ 

plotted  in  Fig.  5.1,  is  2.0  A,  about  a  38%  of  the  discharge  current.  Id  —  5.21  A,  and  a  34% 
of  the  total  ion  production  (about  5.8  A);  Bishaev  and  Kim[l]  estimate,  in  their  experiments, 
relative  current  losses  of  similar  magnitude.  Notice  that  the  plasma  current  lost  in  the  lateral 
walls  implies  a  larger  ion  production  and  eventually  a  larger  discharge  current,  decreasing  thus 
the  thrust  efficiency. 


(5.16) 
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The  electron  energy  deposited  by  conduction  at  lateral  walls,  J  Q^dx,  with 

Qw  Te'^e^we^c 

plotted  in  Fig.  5.1,  is  144  W,  which  represents  a  9%  of  the  electric  power,  ~  1.56  kW;  Ref. 
[1]  estimates  this  energy  loss  around  20%. 

The  structure  of  the  plasma  discharge  in  Fig.  5.1  is  similar  to  the  one  found  when  there  were 
no  lateral  losses.  Fig.  3.2.  It  consists  of  (see  Fig.  2.1)  anode  sheath  (AB),  anode  presheath(BC), 
diffusion  region(CD),  ionization  layer(DH),  acceleration  region(HE)  and  near-plume  (EP).  De¬ 
limiting  points  inside  the  channel  are  xa  —  0,  xc  —  0.5  mm,  X£)  ~  11.5  mm,  Xh  —  14  mm.  The 
regular  sonic  point  is  downstream  the  ionization  layer,  at  xs  —  18.2  mm.  In  the  acceleration 
region,  volume  ionization  and  ion  recombination  at  lateral  walls  are  rather  low,  keeping  the 
ionization  fraction  r]i{x)  almost  constant. 

Notice  in  Eq.(5.14)  for  G  that  the  presence  of  two  terms  for  lateral  losses  gives  more  flexibility 
to  the  balance  necessary  for  a  smooth  sonic  transition  (i.e.  Gs  —  0).  The  drift  of  point  S  away 
from  the  ionization  layer,  when  lateral  losses  are  included,  indicates  that  electron  diffusion  (i.e. 
Joule  heating)  and  lateral  energy  losses  dominate  the  balance  at  the  sonic  transition,  that  is 
we  have 


-U^eTe  ^  ^drrieVl^ 

at  point  S. 

Potential  drops  in  the  different  plasma  regions  are  ^ab  —  1-2  V,  ^cd  —  0.02  V,  cpns  —  —61.9 
V,  (psE  —  —137.6  V,  (pEP  —  —101.5  V,  (j)p^  ~  10  V.  Notice  that  the  potential  drop  in  the 
ionization  layer  is  moderate  and  follows  the  law[37]  4>dh  ~ 

The  thrust,  defined  as 

F  =  {rriinivl  +  Pe)EAc, 

Ss  F  —  90.6  mN  and  the  thrust  efficiency  (using  the  definition  of  Ref.  [37])  is  r/  ~  0.50,  which 
comes  from  the  product  of  utilization  efficiency  ~  0.90,  (electric)  current  efficiency  rjc  —  0.73, 
energy  efficiency  r/g  —  0.97,  and  plume  efficiency  rjp  ~  0.86. 

5.3  Influence  of  control/design  parameters 

Figures  5.3  to  5.6  show  this  influence  for  different  parameters.  These  plots  have  been  obtained 
from  continuation  routines,  departing  from  the  solution  of  Fig.  5.1.  They  can  be  compared 
with  Figs.  3.3  to  3.8  for  zero-losses  at  the  walls.  Sec.  3.6.  Since  wall  losses  are  not  very  large, 
most  features  coincide.  We  observe  that: 

1.  Influence  of  the  mass  flow.  Fig.  5.3:  Discharge  current  and  thrust  increase  linearly  with 
the  mass  flow  -in  agreement  with  experiments  [25,  38,  39,  2]-.  In  addition,  m  is  one  of 
the  few  parameters  affecting  weakly  the  reverse  flow  (i.e.  rjis)-  All  these  features  make 
rh  a  very  suitable  parameter  for  thruster  control.  Although  wall  losses  increase  with  rh 
too,  the  trade-off  for  thrust  efficiency  and  specific  impulse  favors  larger  mass  flows. 
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2.  Influence  of  the  discharge  voltage,  Fig.  5.4:  Both  thrust  and  specific  impulse  grow  like 

-in  agreement  with  experiments  [25,  38,  39,  40]-.  Experiments  [39,  40]  show  that  r] 
presents  a  maximum  for  Va  around  300  V.  The  region  of  positive  slope  for  low  potentials, 
is  not  found  here,  possibly  because  it  corresponds  to  a  regime  of  partial  ionization.  The 
efhciency  decrease  for  large  Va  would  be  due  to  a  larger  acceleration  region  and,  as  a 
consequence,  larger  energy  losses  to  the  walls. 

3.  Influence  of  the  magnetic  field  strength.  Fig.  5.5:  larger  means  an  inhibition  of 
axial  transport,  leading  to  a  larger  ionization  efficiency  and,  eventually,  a  lower  discharge 
current  and  a  larger  thrust  efficiency  favored  by  lower  lateral  losses  too).  Experimentally, 

is  increased  until  a  non-stationary  regime  is  attained[5];  our  results  suggest  that  this 
would  correspond  to  the  limit  rjiB  —  0. 

4.  Influence  of  the  magnetic  field  slope  Fig.  5.6:  The  increase  of  r]  with  L^  indicates  that  it 
is  the  average  magnetic  field  strength  what  influences  basically  the  plasma  response. 

5.  Influence  of  the  cathode  position.  Fig.  5.7:  it  is  relevant  for  short  values  only.  Since  we 
do  not  have  a  clear  criterion  to  choose  xp,  we  suggest  to  avoid  small  values. 

6.  Influence  of  the  channel  length.  Fig.  5.5:  Thrust  and  efficiency  increase  moderately 
with  the  channel  length,  mainly  because  the  plasma  leaves  the  thruster  more  supersonic. 
Experiments  ??  show  the  function  r]{L)  having  a  maximum.  We  could  not  reproduce  it 
here,  possibly  because  wall  losses  have  not  increased  enough  here. 

7.  Influence  of  the  thruster  size.  Fig.  5.9:  Here,  large  wall  losses  make  the  efficiency  to 
present  a  maximum  for  an  intermediate  size.  Notice,  however  that  thruster  operating 
conditions  are  not  optimized  for  each  size. 

Figure  5.10  shows  how  the  pair  of  parameters  (V^,  B^)  should  evolve  to  keep  rjiB  —  —0.05. 
We  expected  that  increasing  both  discharge  voltage  and  magnetic  strength,  rj  would  grow. 
However,  the  increase  of  wall  losses  keeps  the  efficiency  almost  unchanged,  but  thrust 


increases. 
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Figure  5.1:  Axial  structure  of  the  plasma  discharge  for  L  —  25mm,  Lep  —  8.3mm,  A^.  — 
40.1cm^,  Vd  —  300V,  m  —  5.2mg/s,  Tgp  —  2.7eV,  Wap  ~  0.21,  Eyj  —  lOOeV,  i>w{x)  ~  0.17, 
Vnw  —  Vi,  and  n^e/  —  4.3  •  10^®m^^.  (continues  in  Fig.5.2) 
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Figure  5.8:  Influence  of  the  anode  position  (i.e.  the  channel  length  L  —  Xe  —  xa)-  Here,  xa  is 
moved,  while  x^,  Xe  and  Xp  are  kept  fixed.  Rest  of  parameters  as  in  Fig.  5.1. 
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Figure  5.10:  Design  optimization.  Plot  on  the  top-left  shows  the  evolution  of  the  pair  (Va, 
B^)  to  keep  rjis  close  to  0;  rest  of  parameters  as  in  Fig.  5.1.  Rest  of  plots  show  thruster 
performances 


Chapter  6 

EFFECTS  OF  HEAT  CONDUCTION 


6.1  Governing  equations 

Here  we  add  the  effect  of  heat  conduction  into  the  basic  steady-state  model  treated  in  Ch.  3. 
Wall  losses  are  neglected,  i.e.  i>w  —  0.  Again,  for  sake  of  clarity  in  the  reading  we  give  an 
almost  complete  account  of  the  model  equations. 

We  consider  now  the  set  of  equations  We  consider  the  same  stationary  model  of  Chapter  3 
but  including  heat  conduction  effects.  The  governing  equations  are 

d 
dx 

Vn  —  const. 


(AUeV^i)  =  -^{AneV^e)  =  --^{AUnVu)  =  AUel'i 
dx  dx 


Id  „ 

—  —  [AmineV^i)  =  -erie—  +  miUei^iVn, 
dUeTe 

0^  eUe- - - - mel'drieVxe, 

dx  dx 


(6,1) 


Id./  3^  \  ^  dVxe  9  n 

'aIx^  (  -^TeUeVxe  +  Qxe  j  =  +  ^dflefneV^^  -  ViUea^Ei, 


dTe  _  2mel'd 
dx  5neTe 


Qxe- 


6.1.1  Singular  points  and  boundary  conditions 

Solving  Eqs.(6.1)  for  the  spatial  derivatives  we  obtain  again  a  matrix  equation  of  the  form  of 
Eq.(5.11) 


(1  -  M^) 


dY 

dx 


F{Y),  Y  ^  {ne,Vxi,Te,...), 


(6.2) 
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where  now  Y  —  (rig,  Vxi,Te,  q^e, ...)  groups  8  plasma  variables,  F  is  a  regular  function  (different 
from  the  one  in  the  zero  conduction  model),  and  M  is  the  isothermal  Mach  number, 

M  =  ■ 

^/Te/rrii 

For  instance,  the  ion  velocity  follows  Eq.(5.13), 


but  now  function  G  is 


'f^e^xe 


dVxi  ^ 

^1  -  M2’ 


2i'fTlgQxe 
SrigTe  . 


/r.  N  2  din  A 

Vimi[2vxi  -Vn)+  rriiV^i 


(6.3) 


The  eight  boundary  conditions  for  this  model  are  the  seven  ones  of  the  no-conduction  model 
plus  a  condition  for  the  heat  flow  deposited  by  the  electrons  into  the  anode  sheath  (at  point 
B),  already  discussed  in  Sec.  2.6.  Summarizing,  the  eight  boundary  conditions  for  a  solution 
with  supersonic-exit  are 

i)-iv)  m,  VnA  ,  Tep,  and  are  known. 

v)-vii)  At  point  B: 


'^xiB 


=  -\/TeBlmi 


e{^B  -  ^  e^sh  ^  CeB  ^  /  m  \ 

TeB  TeB  V  Fg^  V  / 

QxeB  _  1 

n-eB'^xeBTeB  TeB  2 

viii)  At  point  S  one  has  Ms  —  1  and  Gs  —  0,  that  is 


(6.4) 


Fj^  ^  QxeS-^c 
rUei^dS  2Fjs  STgsFjs 

In  dimensionless  terms,  Eq.  (6.5)  becomes 


ilis) 


mjU* 

rriei^ds 


=  -  r]is  + 


‘^QxeS-^c 

5TesTm 


with  u*  ^—Ris, 
Vn 


(6.5) 


(6.6) 


which  can  be  compared  with  the  corresponding  condition  for  the  non-conduction  problem,  Eq. 
(3.16). 

The  presence  of  an  additional  differential  equation  for  q^e  implies  a  new  mathematical  model, 
with  a  different  sound  velocity  for  the  ion  fluid  and,  therefore,  different  singular/sonic  points. 
The  consequence  is  that  the  model  with  heat  conduction  requires  its  own  integration  algorithms 
and  cannot  be  solved  from  adding  extra  terms  to  the  no-conduction  model  and  performing  a 
parametric  continuation  on  solutions  computed  previously,  nevertheless,  the  integration  pro¬ 
cedure  tries  to  be  similar  to  the  non-conduction  case.  Independent  integrations  are  launched 
from  point  S  towards  points  B  and  P.  As  before,  region  BS  poses  the  main  restrictions  to  the 
existence  of  solutions.  An  additional  difficulty  now  is  that  the  solution  launched  from  S  must 
meet  condition  (6.4)  on  qxeB,  which  requires  iteration  on  the  new  ’initial’  parameter  qxes-  In¬ 
deed,  due  to  convergence  problems,  it  is  easier  to  launch  solutions  from  points  B  and  S  and  to 
match  them  in  an  intermediate  point. 
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6.1.2  Restriction  on  the  ion  backflow 


One  new  feature  of  the  conduction  model  is  the  restriction  of  the  relative  ion  back- flow,  rjis, 
to  a  narrow  range.  The  matching  of  the  quasineutral  solution  with  the  anode  sheath  requires 
(i)  the  anode  sheath  to  be  ion-attracting  and  (ii)  the  ion  back-flow  to  accelerate  towards  the 
anode.  The  first  condition  means  >  0.  From  Eq.(3.12),  the  second  condition  implies  that 
Gb  >  0  in  Eq.(3.12),  or,  using  Eq.(6.3)  with  Ui  ~  0, 


‘^Qx 


dnpTpVj 


<  1. 


(6.7) 


Using  now  Eqs.  (3.11)  and  (2.61),  the  two  conditions  yield  the  double  restriction 

0<ln-^^<3,  (6.8) 

4\VxeB\ 

or,  in  terms  of  the  ion  back  current, 

5 . 10-3  <  <  0.114,  (6.9) 

Id 

with  Ii  —  eniVxiAp.  However,  difficulties  found  in  the  numerical  convergence  around  point  B 
near  the  limit  (^ab  —  0,  have  forced  us  to  exchange  that  limit  for  the  close  one,  q^eB  —  0,  where 
e4^AB/TeB  —  0.5.  This  moves  the  lower  bound  of  the  ion  back-current  to  llwl/^d  ^  8  •  10“^. 


6.2  Solution  for  zero- ionization  regions 


Since  the  ionization  region  is  relatively  thin,  it  is  of  interest  to  know  how  is  the  solution  in  the 
regions  where  ionization  is  negligible.  Making  r',  =  0  in  Eqs.  (6.1)  one  has 


AuaVa  —  const  (cK  =  i,  e,  n), 
1  2 

-rriiV^i  +  —  const, 

—  e$  H — =  const, 

2 


dVxi 

rriineVxi— - h 

ax 


dnpTp 

dx 


-vampnpVxe- 


(6.10) 


These  equations  together  with  the  equation  for  dTp/dx  yield  a  second  order  differential  equation 
for  Tp  and  q^p  =  qxe/npVxpTp\ 

1  dTp  _  ^ 

Tp  d^  ~ 

(  _  1  ')  ^dxe  _  1  I  f  ^  -2  ,  ^  dA 

\mivT  )  d^  Vs  rriivlJ^'^''  5v  rriivlJ^''''  A  d^' 

dx  Tp 

d^  VdmpVxp ' 


(6.11) 
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The  (dimensionless)  auxiliary  variable  ^  is  used  for  convenience  and  to  show  the  spatial  scale 
of  the  problem. 

The  equation  for  dq^e/d^  shows  that  a  regular  sonic  transition  is  possible  now  even  without 
ionization.  Sonic  point  S  is  regular  if  Gs  —  0  becomes  here 

5  5 

feS  =  2’  that  is  QxeS  ^  ^{neVxeTe)\s,  (6.12) 

a  result  that  can  be  obtained  from  Eq.(6.5)  also.  Therefore,  now,  a  regular  sonic  point  can  exist 
even  in  a  region  of  negligible  ionization.  This  is  a  novelty  with  respect  to  the  non-conduction 
model. 

A  particular  solution  of  these  equations  for  A  —  const  is 

Te,  $  =  const,  Ua,,,  =  0,  (6.13) 

which  corresponds  to  the  asymptotic  behavior  of  the  diffusion  region,  point  C. 

6.3  Plasma  structure  and  thruster  performances 

Figure  4  shows  the  axial  profiles  of  main  plasma  variables  for  an  SPT-100  type  of  thruster. 
Design  and  control  parameters  are  similar  to  the  case  of  Fig.  3.  The  plasma  structure  presents 
similar  trends  to  the  no-conduction  solution  of  Ref. [37]  (with  no  lateral  losses),  but  heat  con¬ 
duction  smoothes  the  temperature  profile  in  two  ways:  first,  the  peak  temperature  is  reduced 
from  about  90  eV  to  65  eV;  second  there  is  a  smoother  transition  between  the  diffusion  and 
ionization  regions,  leading  to  a  shorter  diffusion  region,  and  a  larger  ionization  layer.  The  re¬ 
duction  in  the  maximum  temperature  is  not  enough  to  agree  with  experimental  values  where, 
for  Vd  ~  300  V,  maximum  temperatures  are  in  the  range  20  —  30eF.  Therefore,  wall  losses  need 
to  be  considered  still. 

Plasma  equations  and  boundary  conditions  seem  to  force  an  almost  monotonic  profile  of  q^e 
between  B  and  S.  Then,  the  heat  flux  at  point  S,  qxes^  is  large.  In  the  present  model,  Fq.  (6.3) 
shows  that  a  regular  sonic  transition  implies  to  balance  Joule  heating,  conductive  heat  flow, 
and  ionization  losses.  The  large  value  of  q^es  means  this  balance  to  take  place  mainly  between 
heat  conduction  and  Joule  heating, 

5  ^ 

Qxe  '^'f^eTe'Vxei 

and,  therefore,  point  S  is  outside  (downstream,  in  fact)  the  ionization  layer.  In  turn,  the  point 
of  maximum  temperature,  where  qxe  —  0,  is  going  to  be  downstream  of  point  S,  what  justifies 
that  we  could  not  find  valid  solutions  with  the  point  of  maximum  temperature  well  inside  the 
channel. 

For  the  case  of  Fig.  4,  Fig.  5(a)-(c)  measure  different  phenomena  in  the  electron  dynamics. 
This  is  aimed  to  evaluate  the  validity  of  the  hypotheses  supporting  our  model.  Figure  5(a) 
shows  the  different  contributions  to  the  total  collision  frequency  of  electrons,  Fq.(2.57).  Bohm 
diffusion  provides  the  main  contribution  from  the  ionization  layer  towards  the  exit,  whereas 
e-i  collisions  (instead  of  the  e-n  collisions)  dominate  in  the  diffusion  region.  This  last  feature 
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is  justified  by  the  high  plasma  density  and  low  temperature  there,  and  needs  further  investi¬ 
gation.  Figure  5(b)  shows  the  axial  profiles  of  the  ionization  frequency,  with  its  peak  marking 
the  ionization  region,  and  the  effective  axial  frequency  for  electron  diffusion,  Eq.(134);  the 
dominance  of  Bohm  diffusion  makes  oc  B{x),  approximately,  in  most  of  the  channel.  Figure 
5(c)  shows  that  the  electron  mean  azimuthal  energy  (main  contribution  to  the  electron  mean 
kinetic  energy)  is  totally  negligible  compared  to  the  thermal  energy,  except  perhaps  within  a 
small  region  near  the  anode. 

Performance  parameters  for  the  case  plotted  in  Fig.  4  are:  F  ~  90.6mN,  ~  4.51A, 

{id  —  1.17),  and  rj  ~  57.5%.  The  influence  of  the  plasma  discharge  Vd  and  the  magnetic  field 

on  thruster  performances  have  been  computed  from  parametric  continuation  of  the  solution 
of  Fig.  4.  The  consequence  of  the  narrow  interval  for  rjis  in  Eq.(6.9),  is  a  narrow  range  too 
of  the  control  parameters  where  our  model  yields  stationary  solutions  (with  large  ionization). 
As  an  example.  Fig.  6(a)  shows  the  band  of  of  parameters  {Vd,  B^),  where  solutions  are  found 
is  dramatically  small;  only  variations  of  about  10%  are  permitted.  Figures  6(b)-6(d)  plot  the 
corresponding  performance  bands  for  thrust,  discharge  current,  and  efficiency;  for  F,  the  band 
reduces  practically  to  a  line.  These  plots  show  that  larger  thrust  and  efficiency  are  obtained 
with  larger  Vd  and  B^. 
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Figure  6.1:  Model  with  heat  conduction  (and  no  wall  losses).  Axial  structure  of  the  plasma 
discharge  for:  L  —  25mm,  Lep  —  21.1mm,  Ac  —  45cm^,  —  300V,  m  —  5.2mg/s,  Tgp  —  2.5eV, 

IGap  ~  0.067. 
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Figure  6.2:  Model  with  heat  conduction  (and  no  wall  losses).  Same  case  of  Fig.  5.  (a) 

Contributions  of  different  processes  to  the  electron  collision  frequency  r'g-  (b)  Axial  diffusion 
frequency  (1/^  =  and  ionization  frequency  along  the  channel,  (c)  Ratio  of  electron 

mean  kinetic  energy  to  electron  temperature. 
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Figure  6.3:  Model  with  heat  conduction  (and  no  wall  losses),  (a)  Band  of  input  parameters 
{Vd,  B^)  where  stationary  solutions  have  been  found;  rest  of  parameters  as  in  Fig.  4.  Line  ’a’ 
corresponds  to  Gb  —  0,  and  line  ’b’  to  (pAs  —  0.  Bands  for  (b)  thrust,  (c)  discharge  current, 
and  (d)  efhciency,  corresponding  to  the  (14,  B^)  band. 


Chapter  7 

SELF-OSCILLATIONS 


7.1  Perturbation  equations 

The  oscillatory  behavior  of  Hall  thrusters  is  here  approached  by  investigating  the  stability  of 
the  steady-state  plasma  structure  studied  in  previous  chapters.  The  linear  stability  analysis 
is  based  in  the  plasma  response  to  small  temporal  perturbations.  Equations  and  results  are 
presented  here  for  the  no-conduction  model  with  wall  losses  of  Chap.  5. 

Since  temporal  derivatives  affect  to  the  perturbation  terms  only,  it  is  convenient  to  solve 
plasma  equations  (2.55)  for  the  spatial  derivatives.  This  yields 
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with 


P^T 
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[2  3 

G  =  VxgG  Uj  —CXiEf  “h  Tg  '^fTT'i'Vxeip'^xi  ^n) 


^dX^e^xe  \ 


+  -T^weTg  +  -rUivliV 
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+  ’zfTli'Vxg  I  Vxi' 
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5  dt 


- 


dt 


The  equation  for  the  area  variation  in  the  plume  must  be  added: 


da 

‘dx 


f  0,  X  <  L^ 

2cp/d,  L  <  X, 


(7,2) 


(7.3) 


with  a  —  \nd. 


7.2  Model  of  linear  perturbations 


To  be  definite,  let  us  assume  a  small  perturbation  in  the  channel  discharge  voltage, 

Vd{t)^Vdo  +  u{Vd^e-^^^], 

with  Vdi  <C  Vdo  and  uj  —  Urg  +  icoim  the  (complex)  frequency  of  the  perturbation  mode.  Then, 
in  the  linear  approximation,  the  plasma  response  is 

Y{x,t)  —  Yoix)  +  3fi|Ti(a;;  w)  exp(— ia;t)|, 

with  Y{x,t)  the  vector  of  plasma  variables,  Yq{x)  the  solution  obtained  from  the  stationary 
model  of  Chapter  5  (subscript  0  is  added  now  to  stationary  variables),  and  Ti(a;;a;)  the  linear 
perturbation  response,  which  is  proportional  to  Vdi- 
From  the  expansion  of  Eqs.(7.1)  and  the  substitution 

d/dt  —in;. 


the  linear  equations  for  each  perturbation  mode  in  oj  are 
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(7.4) 

(7.5) 

(7.6) 

(7.7) 

(7.8) 

(7.9) 
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dvni  dv„o 
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plus  the  equation  for  the  area  variation 


X  <  L, 


-  {aiVxio  A  Vxii)^  +  2^  L  <x 
ox  ao 


(7.15) 


(7.16) 


with  gq  —  In  do,  =  doo^i-  Finally,  the  linear  expansions  for  Vdi,  r'ji,  i^wi,  and  r'.^ei  niust  be 
added  to  these  equations. 


7.2.1  Boundary  conditions 

Like  the  stationary  equations,  the  perturbation  equations  are  singular  at  Pq  —  0.  The  appro¬ 
priate  set  of  seven  boundary  conditions  for  the  perturbation  problem  is: 
i)-iv)  thiA,  VniA,  Teip,  and  <j)ip  =  -  V^i  are  known. 

v)  To  fulfill  hgiB  ^  heQBi  the  singularity  at  point  B  for  the  perturbed  response  cannot  be 
of  higher  order  than  for  the  stationary  response.  Therefore,  it  must  be 


PlB  =  0. 

vi)  Taking  (ji^A  =  0,  from  Eq.(3.11),  the  perturbation  of  the  sheath  potential  is 


^(t>lB  —  ^4^0B7^ - H  TgQB  I 

eOB 


■  iOB  eOB  ' 


(7.17) 


(7,18) 


vii)  For  the  perturbed  solution  to  be  regular  at  point  S,  the  right  side  of  Eq.(7.4)  must  be 
zero.  Considering  that  Pqs  =  Gqs  =  0,  the  regularity  condition  at  point  S  is 

0  =  [Gi  -  =  G^s  -  (7-19) 

L  ox  is  OPo/Ox\s 
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7.2.2  Integration  procedure 

The  perturbation  of  point  S,  xsi,  is  part  of  the  solution.  From  the  expansions  of  =  0  and 
Gs  =  0,  one  has 

-Pis  Gis  ,  . 

xsi  =  -op  /o  I  =  I  •  (7-20) 

dPo/dx\s  dGo/dxls 

The  compatibility  of  these  two  relations  is  assured  by  condition  (7.19). 

The  perturbation  of  the  discharge  current,  r,^!,  is  part  of  the  solution  also.  It  satisfies 

rji(a;)  -  Tei{x)  =  Tdi  =  const,  (7-21) 

However,  the  mass  flow  perturbation  changes  along  the  channel, 

Tj^a;)  +  r„i(a;)  =  r„i(a;)  const. 

The  numerical  integration  of  this  perturbation  problem  is  performed  again  with  a  Runge- 
Kutta  algorithm  from  point  S  towards  points  B  and  E.  However,  there  appear  new,  delicate 
issues  in  the  region  BS,  related  to  singular  points  B  and  S.  These  are 

1)  To  obtain  the  spatial  derivatives  of  perturbed  variables  at  point  S  requires  to  know  the 
second  derivatives  of  Gq  and  Pq. 

2)  These  second  derivatives  must  be  used  now  on  the  0th  order  solution  as  well,  for  a  smooth 
matching  of  the  spatial  profiles  in  the  vicinity  of  point  S. 

3)  A  perturbation  solution  launched  from  point  S  is  unbounded  (i.e.  tends  to  infinity)  at 
point  B.  Therefore,  it  must  be  stopped  before  point  B  and  be  matched  to  an  independent 
solution  launched  from  point  B,  at  an  intermediate  point  (here  chosen  as  Mq  =  —0.2).  This 
increases  the  number  of  ’fundamental  modes’  (each  one  proportional  to  a  ’fundamental  param¬ 
eter’)  constituting  the  general  solution,  but,  since  the  problem  is  linear,  it  avoids  completely 
any  subsequent  iteration. 

The  solution  of  the  perturbation  problem  is  thus  written  in  the  form 

Y,{x)  =  +  Y^buYnix),  (7,22) 

j  I 

where 


Gi  —  [Cll ,  . . . ,  C19]  l^elS  T  ^nlS T  ^xelS T  ^nlS  T  ^elS T  4^1S T  ^elB T  ^xelB ]  (7.23) 

is  a  vector  of  unknown  parameters,  [611,612]  =  [vniAi^miA  known  parameters,  and  Y^j  and 
Y^i  constitute  eleven  ’fundamental  modes’,  which  depends  on  a)  and  the  vector  Cq  of  parameters 
defining  the  stationary  solution  (tildes  mean  dimensionless  magnitudes). 

Next  we  impose  the  continuity  of  the  solution  at  the  intermediate  point  and  the  values  Tgip 
and  (j)ip  for  the  temperature  and  the  potential  at  point  P.  In  total,  these  constitute  nine  linear 
relations  that  we  can  write  in  the  form 


MCi  =  Hi,  (7.24) 

where:  Hi  includes  five  zeros  and  the  perturbations  of  four  control  parameters:  VniA,  ^miA-, 
Tgip  and  ^ip;  and  M  is  a  square  matrix,  obtained  from  the  fundamental  modes,  which  depends 
on  a)  and  Cq. 
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7.3  Self-excited  modes 

Linear  self-excited  modes  are  singular  solutions  of  the  perturbation  problem  for  homogeneous 
boundary  conditions,  i.e.  Bi  =  0.  The  frequencies  of  the  self-excited  modes  are  the  solutions 
of  the  global  dispersion  relation 

D(uj;  Co)  =  det  M  =  0.  (7-25) 

The  self-excited  modes  are  unstable  if  >  0. 

Fig.  7.1  shows  the  frequencies  of  the  main  five  self-excited  modes  corresponding  to  the 
steady-state  solution  of  Fig.  5.1.  Notice  that  only  Mode  1  is  unstable  for  this  case.  Self-excited 
modes  are  difficult  to  obtain  since  the  zeros  of  detD  are  very  sharp.  This  makes  advisable  to 
solve  both  Dre  =  0  and  =  0  instead  of  just  detD  =  0;  indeed,  Noguchi[41]  failed  to  find 
unstable  modes  because  he  used  this  last  method. 

The  five  self- frequencies  of  Fig.  7.1  are: 

-  Mode  0:  a;(kHz)  =  0  —  7.2i7.2 

-  Mode  1:  a;(kHz)  =  45.2  +  17.2i 

-  Mode  2:  a;(kHz)  =  142.3  -  23.2i 

-  Mode  3:  a;(kHz)  =  51.7  —  75. 3i 

-  Mode  4:  a;(kHz)  =  119.5  —  90. 4i 

Figures  7.2-7. 8  show  the  spatial  profiles  of  some  perturbation  variables  for  the  five  modes, 
and  contour  plots  for  Modes  1  and  2.  These  two  modes  seem  to  be  the  main  ones  with  respect  to 
stability  (Mode  2  has  been  found  unstable  for  other  steady-state  solutions).  The  five  self-excited 
modes  present  several  similarities.  First,  in  all  modes,  unstable  and  stable,  perturbations  are 
stronger  in  the  ionization  layer.  This  is  the  consequence  of  spatial  gradients  being  larger  there. 
Large  ionization  fluctuations  are  found  also  in  the  dynamic  response  to  regular  perturbations, 
-that  is  solutions  with  Bi  ^  O  in  Eq.  (7.24)-.  Therefore,  ionization  waves  should  not  be 
considered  the  mark  characterizing  self-excited  modes,  as  suggested  by  some  people.  Other 
common  features  of  self-excited  modes  are  that  ion  and  neutrals  acoustic  waves  are  launched 
in  both  directions  within  the  subsonic  region  of  the  channel,  whereas  the  acceleration  region 
and  the  plume  seem  to  have  a  minor  role  in  the  dynamic  response. 

Contour  plots  in  Figs.  7.4  and  7.6  show  that  neutral  perturbations  are  forward  travel¬ 
ing  waves  (indeed  they  are  standing  waves  mounted  on  the  stationary  neutral  flow),  while 
ion  perturbations  form  a  backward  traveling  wave.  Electron  perturbations  are  controlled  by 
quasineutrality  and  fulfill  Eq.  (7.21).  The  main  differences  between  Modes  1  and  2  are  observed 
in  the  wavelengths  and  phase  velocities  of  the  ion  and  neutral  waves.  This  seems  to  indicate 
that  the  self-excitation  is  related  to  the  coupling  between  ionization  and  acoustic  waves  in  the 
diffusion  region,  with  a  secondary  role  for  the  acceleration  region.  However,  further  work  is 
necessary  to  understand  the  sequence  of  mechanisms  that  feed  the  self-excited  modes. 

Figures  7.9-7.11  show  the  influence  of  Va  and  rh  in  the  self- frequency  of  Mode  1.  Relative 
variations  are  of  order  one  and  the  instability  persists  in  the  parametric  range  shown  there. 
No  general  conclusions  should  be  drawn,  from  this  particular  case,  about  the  stabilizing  or 
unstabilizing  effect  of  those  parameters. 

Finally,  the  comparison  of  self-excited  modes  of  steady-state  solutions  with  and  without 
wall  losses  does  not  show  significant  differences. 
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Figure  7.1:  First  self-excited  modes  for  the  stationary  solution  of  Fig.  3.2.  Solid  lines:  D^e  =  0; 
dash-dot  lines:  =  0.  Modal  frequencies  are  given  by  the  intersection  of  both  lines  Dashed 

lines  are  contour  lines  for  ID  I  constant  and  small. 


x(mm} 


x{mm) 


Figure  7.2:  Mode  0  (of  Fig.  7.1):  Spatial  profiles  of  different  perturbation  magnitudes  [say 
'Ki(a;,  t)  =  F^i(a;)  exp(— iwit)]  normalized  with  T^x/^do  =  1-  Solid  line:  envelope  of  perturbation 
amplitude  [i.e.  |'Ki(a;)|];  Dashed  line:  perturbation  at  a  fixed  time  [i.e.  (^)}|] 
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Figure  7.9:  Evolution  of  the  complex  frequency  of  Mode  1  with  The  rest  of  0th  order 
parameters  as  in  Fig.  3.2. 


Figure  7.10:  Evolution  of  the  complex  frequency  of  Mode  1  with  V^.  The  rest  of  0th  order 
parameters  as  in  Fig.  3.2. 


Figure  7.11:  Evolution  of  the  complex  frequency  of  Mode  1  with  m.  The  rest  of  0th  order 
parameters  as  in  Fig.  3.2. 


Chapter  8 
SUMMARY 


This  work  has  presented,  through  successive  stages  grouped  in  Chaps.  2  to  7,  a  semi-analytical 
macroscopic  model  of  the  plasma  discharge  in  a  Hall  thruster.  In  Chap.  2  we  have  formulated 
the  general  model  for  the  plasma  discharge  both  inside  the  thruster  and  in  the  near-plume. 
Plasma  equations,  collisional  processes,  simplifying  hypotheses,  spatial  scales,  and  other  issues 
have  been  discussed  and  evaluated. 

Chapter  3  presents  what  can  be  called  the  ’core  model’,  which  ignores  electron  heat  con¬ 
duction  and  the  interaction  of  the  plasma  with  the  channel  lateral  walls.  Since  both  thruster 
channel,  near-plume,  and  neutralizer  are  included  in  the  model,  we  could  analyze  the  complete 
response  of  the  plasma  in  terms  of  actual  control/design  parameters.  The  main  achievement 
of  the  core  model  is  the  reproduction  of  the  different  regions  of  the  rich  plasma  structure  ob¬ 
served  experimentally,  and  the  understanding  of  many  of  the  discharge  basic  physics.  Other 
accomplishments  here  have  been  (i)  to  show  that  optimum  operation  conditions  correspond  to 
solutions  with  an  internal  sonic  transition,  and  (ii)  to  find  that  the  domain  where  stationary 
solutions  exist  is  limited. 

In  Chap.  4  we  discussed  a  radial  model  for  the  interaction  of  the  plasma  with  the  dielectric 
walls  of  the  thruster  channel.  The  model  includes,  in  a  consistent  way,  wall-emitted,  plasma- 
trapped  electrons,  explains  the  transition  to  a  regime  with  charge-saturated  sheaths,  and  solves 
the  presheath/sheath  radial  structure  totally.  Finally,  losses  of  plasma  current  and  energy,  and 
gains  in  electron  azimuthal  momentum  (leading  to  near- wall  conductivity),  in  the  channel,  have 
been  evaluated. 

Chapters  5  and  6  were  dedicated  to  improve  the  core  model.  In  Chap.  5  we  included  the 
effects  of  the  radial  plasma/wall  interaction  into  the  axial  model.  Results  obtained  until  now 
show  very  large  losses  of  plasma  current  and  energy,  unless  some  mechanism,  not  well  identified 
yet,  inhibits  the  radial  fluxes  to  the  walls.  On  the  contrary,  near-wall  conductivity  is  found 
marginal,  for  inhibited  flows,  at  least.  However,  few  cases  have  been  run  still  and  further 
research  is  needed  before  reliable  conclusions  can  be  formulated. 

In  Chap.  6  we  modified  the  core  model  by  considering  heat  conduction  effects  (but  radial  in¬ 
teraction  was  ignored).  The  differences  in  the  mathematical  formulation  for  the  heat-conduction 
model  has  required  to  adapt  the  integration  procedure.  The  new  stationary  solutions  present 
smoother  temperature  profiles,  more  according  to  experimental  ones,  but  heat  conduction  is 
found  insufficient  by  itself  to  reduce  the  peak  temperature  to  experimental  values.  The  general 
plasma  structure  is  maintained  with  a  shorter  diffusion  region  and  the  sonic  transition  away 
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from  the  ionization  layer  outer  bound.  Finally,  the  most  worrying  conclusion  of  this  conductive 
model  is  the  dramatic  reduction  in  the  parametric  domain  of  stationary  solutions  (of  known 
ones,  at  least). 

In  Chap.  7  we  undertook  a  linear  stability  analysis  of  the  steady-state  plasma  structure. 
The  analysis  is  classical  but  the  formulation  and  solving  of  the  perturbation  model  is  rather 
sophisticated  due  to  the  rich  stationary  structure  and  the  presence  of  singular  points.  Lon¬ 
gitudinal  self-excited  modes,  some  of  them  are  stable  and  other  ones  unstable,  are  found  in 
the  same  frequency  range  than  experiments  (10-100  kHz),  Self-oscillations  are  far  from  being 
simple  and  well  understood.  They  consist  of  local  ionization  waves  and  ion  and  neutral  acoustic 
waves  counter-traveling  along  the  subsonic  part  of  the  channel.  Coupling  of  the  different  funda¬ 
mental  waves  takes  place  at  the  ionization  layer  and  the  anode  sheath,  and  through  collisional 
processes.  Electron  perturbations  are  governed  by  quasineutrality,  and  the  regions  of  super¬ 
sonic  acceleration,  including  the  plume,  have  a  marginal  role  in  these  modes.  An  unfinished 
parametric  investigation  have  not  shown  ways  to  reduce  effectively  the  unstable  modes  yet. 

Different  parts  of  this  work  have  been  presented  in  recent  International  Conferences  and 
published  in  technical  papers  and  in  Physics  of  Plasmas  [42,  37,  43,  36,  44]. 

Future  research  should  deal  with: 

1)  A  model  of  radial  interaction  with  partial-trapping  of  secondary  emission  and  cylindrical 
effects.  Partial-trapping  is  expected  to  obstruct  charge-saturation  of  the  radial  sheath  and, 
therefore,  to  reduce  energy  losses.  Cylindrical  effects  should  produce  an  asymmetry  between 
the  two  thruster  walls  with  possible  effects  on  electron  trapping  and  flow  inhibition. 

2)  A  unified  model  with  heat  conduction  and  losses  at  lateral  walls.  This  means  to  add  the 
radial  model  of  Chap.  4  (or  an  improved  version)  into  the  conductive  model  of  Chap.  6. 

3)  A  comparative  study  of  the  results  of  these  ID  macroscopic  models  with  those  of  2D 
PIC-based  models.  This  comparison  could  provide  relevant  clues  to  evaluate  and  improve  both 
types  of  models  and  to  use  them  in  combination  for  thruster  design  and  analysis. 

4)  Progresses  in  the  study  of  thruster  oscillations.  This  should  include  both  a  better  under¬ 
standing  of  the  unstable  modes,  and  investigation  of  whether  the  parametric  domain  with  no 
stationary  solutions  corresponds  to  thruster  oscillatory  behaviors  observed  experimentally. 


Appendix  A 

Expansion  around  sonic  points 


A.l  Stationary  model 

A.  1.1  Singular  sonic  point 

At  a  singular  sonic  point  one  has  Pq  =  0  and  Gq  ^  0.  Around  this  point  and  for  djdt  =  0, 
Eqs.(7.1)  yield  the  expansion 


GqcIx  =  PodhgQ, 

dhgQ  —  dv  xio  /  V  xio  j 

dTgO  —  {Tgo/Uxio  nT'i'Vxio}dVxiO- 

Taking  Vxio  as  independent  variable,  and  differentiating  once  more,  one  has 


d'^x  ~ 

5Go 

and  the  velocity  profile  around  the  singular  point  is 


2 

xiO') 


^Vxio  = 


(A.l) 


(A.2) 


The  profiles  of  Tgo  and  hgo  are  obtained  from  the  above  linear  relations  with  Vxio-  From  Eq.(A.2), 
a  valid  spatial  solution  with  points  B  and  S  singular,  must  have  <  0  and  Gs  >  0. 


A.  1.2  Regular  sonic  point 

In  a  regular  point  S,  one  has  Gqs  =  0  and  Pqs  =  0,  so  the  derivatives  of  the  plasma  variables  in 
Eqs.(7.1)  are  not  determined  and  the  numerical  integration  cannot  proceed  directly  from  point 
S.  To  determine  them,  it  is  convenient  to  use  the  auxiliary  variable  ^  defined  by 

Xdx/d^  =  Xx'  =  Po'^a^eo,  (A. 3) 

with  A  any  constant.  Then,  the  equation  for  in  Eqs.(7.1)  becomes  (for  A{x)  =  const  around 
point  S) 

=  ^0  (A. 4) 
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and  Eqs.(7.1)  can  be  written  as 


=  -VxioKo  +  (A.5) 

TeO  =  -rriiVxiov'^iO  -  -  Vxio)  -  T^dorrieV^eo^  x',  (A.6) 

=  -VxeoKo  +  (A.7) 

=  T^o  +  TeoKo  +  mei^doVxeox',  (A.8) 

<0  =  0,  (A.9) 

Vnon'r^o  =  -neoViQx'.  (A.IO) 


Notice  that  point  S  is  now  a  ’stationary  point’  of  the  above  equations  in 
The  expansion  of  the  RHS  of  Eqs.  (A. 3)  and  (A. 4)  around  point  S  is 

-  Ko  -  ^mvxiov'^io  +  HOT, 

3  2 

Afc'u  ~  -  2l'jo)»ieO<'i|o  +  (I'iO  -  I'«.0  +  -  mevl^v'^ 


(A.ll) 


3  2  1  /3  \  2 

T  '^'^i'^xeoiTnO  ‘^'^xio)  T  T  ^eO  ^jo  xiO^ xeO  ^eO  1  ^wO  '^'^eO^weO  HOT, 

(A.12) 


or,  in  matrix  notation. 


~  Ml 


+  HOT. 


(A.13) 


The  derivatives  of  the  diverse  frequencies  can  be  written,  in  matrix  form,  as 


'do  _ 


= 


(A.14) 


Then,  Eqs.  (A.5)-(A.6)  provide  a  matrix  relation 


Neglecting  now  HOT  in  Eq.  (A.13),  one  has,  at  point  S, 


Ko\^ 


(A.15) 


(A.16) 


{il^i  +  +  M2^3- 


(A.17) 
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The  two  possible  slopes  of  dheo/x'\s  correspond  to  the  eigenvectors  of  matrix  A.  These  eigen¬ 
vectors  are  real  for 


(a  —  -I-  46c  >  0, 

so  only  the  parametric  domain  verifying  this  condition  must  be  studied. 


(A.18) 


A. 2  Perturbation  model 

A. 2.1  Singular  sonic  point 

To  depart  from  point  B,  the  quantity 


—  Po'^xeO 


dYi 

dx 


E,B 


must  be  known,  which  means  to  solve  the  indetermination  Pi/Pq  in  Eq.(7.4).  Expanding  the 
perturbation  equations,  the  problem  reduces  to  solve  three  linear  equations  for 


—  Q  _  ( ^el  q 

^a:*l  — 

^el  “  ~  ~  —  Tgo/igj, 


(A.19) 

(A.20) 

(A.21) 


with 


'^xioOoi 


T' 

-‘-eO 


i'^liO  ~  Teo)Go- 


(A.22) 


A. 2. 2  Regular  sonic  point 

Equations  need  to  be  expanded  twice  to  determine 

Y;  =  dYi/dx\s. 

Differentiating  in  Eq.  (5.5)  for  hgo,  one  has 

PoVxeoKo  +  Po'^xeoKo  +  Pov'xeoKo^  =  ^0 


and,  particularizing  at  point  S, 

Pos'^xeoshgQs  =  Gq^. 

Differentiating  in  Eq.(A.23)  and  particularizing  at  point  S,  one  has 

‘^Pos'^xeOshf.Qs  =  —  {Pos'^xeOS  +  ^Pos'^xeOs)KoS 


(A.23) 


(A.24) 


(A.25) 
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Finally,  differentiating  Eq.(7.4), 

Po'^xeohgi  +  PoVxeohgi  +  -Po'^a:e0^el  “ 

G'l  —  {PoVxel  +  PlVxeo)Ko  ~  {Po'^xel  +  -Pl^^a^eO  +  Po^xel  +  Pl'^xeo)Ko 


and  particularizing  at  point  S,  one  has 

Pos'^xeoshf.is  =  Gig  —  PlsVxeOshf,Qg  —  {Pos'^xelS  +  Pls'^xeOS  +  PlsVxeOs)^eOS-  (A. 26) 

Equations  for  Pq,  Pq  ,  Gq,  Gq,  P{,  and  G'l  include  first  and  second  derivatives  of  the  different 
plasma  variables.  Once  all  the  equations  are  added  we  have  three  set  of  linear  equations  to 
determine  Y^g,  Y^'g,  and  Y{g. 
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